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Abstract 

The heat kernel expansion is a very convenient tool for studying one-loop diver- 
gences, anomalies and various asymptotics of the effective action. The aim of this 
report is to collect useful information on the heat kernel coefficients scattered in 
mathematical and physical literature. We present explicit expressions for these co- 
efficients on manifolds with and without boundaries, subject to local and non-local 
boundary conditions, in the presence of various types of singularities (e.g., domain 
walls) . In each case the heat kernel coefficients are given in terms of several geometric 
invariants. These invariants are derived for scalar and spinor theories with various 
interactions, Yang-Mills fields, gravity, and open bosonic strings. We discuss the re- 
lations between the heat kernel coefficients and quantum anomalies, corresponding 
anomalous actions, and covariant perturbation expansions of the effective action 
(both "low-" and "high-energy" ones). 
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1 Introduction 



It was noted by Fock [207] in 1937 that it is convenient to represent 
Green functions as integrals over an auxiliary coordinate (the so-called "proper 
time" ) of a kernel satisfying the heat equation. Later on Schwinger [379] recog- 
nised that this representations makes many issues related to renormalization 
and gauge invariance in external fields more transparent. These two works 
introduced the heat kernel to quantum theory. DcWitt made the heat kernel 
one of the main tools of his manifestly covariant approach [147,148,149,150] to 
quantum field theory and quantum gravity which became dominant for many 
years. 

Heat kernel is a classical subject in mathematics [271] ^ . Asymptotics of 
the heat kernel are closely related to the eigenvalue asymptotics found by 
H. Weyl [418,419] and studied further in [325,326]. The problem, as it was 
formulated by Kac [286], reads: "Can one hear the shape of a drum?". In 
other words, this is the problem of recovering geometry of a manifold from 
the spectrum of a natural differential operator. Heat kernel coefficients proved 
very useful in this context. On the other hand, the heat kernel is also an 
adequate tool to study the index theorem of Atiyah and Singer [22,236, 18]. 

By about 1990 the heat kernel expansion on manifolds without bound- 
aries or with boundaries and simplest local boundary conditions on them was 
well understood. Also, the heat kernel became a standard tool in calcula- 
tions of the vacuum polarisation, the Casimir effect, and in study of quantum 
anomalies. Later on, progress in theoretical physics, especially in string the- 
ory and related areas, and parallel developments in mathematics made this 
field highly specialised. New results on non-standard boundary conditions (as, 
e.g., containing tangential derivatives on the boundary or non- localities) or on 
non-standard geometries (domain walls) were scattered in large amounts of 
physical and mathematical literature. The aim of this report is to present a 
unifying approach to the heat kernel expansion and to supply the reader with 
a "user-friendly" guide to the field. The main idea which we shall pursue is 
the universality of the heat kernel. A single calculation (though, sometimes, 
quite involved) may help in a large variety of applications regardless of such 
details as spin, gauge group, etc. As well, just a single universal object in fact 
describes counterterms, anomalies, some asymptotics of the effective action, 
and much more. 

To illustrate the use of the heat kernel in quantum field theory let us 
consider the generating functional for the Green functions of the field (j) in the 
path integral representation 



A historical survey of the mathematical Uterature on the heat kernel expansion 
can be found in [245]. 
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The heat kernel methods arc almost exclusively used for the one-loop calcula- 
tions. In this approximation it is enough to expand the action L to up to the 
quadratic order in quantum fluctuations 0. 

£ = £ei + (0,J) + (0,i?0), (1.2) 

where £ci is the action on a classical background, (., .) denotes an inner product 
on the space of quantum fields. Usually, this inner product is just an integral 
over the underlying space. For real one-component scalar fields it reads 

(01, 02) = / rX^Mx)cl>2{x) . (1.3) 

The linear term in (1.2) contains in general contributions from the external 
sources of the field 4> and from the first variation of the classical action. If 
the background satisfies classical equation of motion for the field 0, the latter 
part of the linear term vanishes, though the former one (external sources) 
should be kept arbitrary if one wishes to study correlation functions of 0. We 
stress that the background and quantum fields may be of completely different 
nature. For example, it is a meaningful problem to consider pure quantum 
scalar fields on the background of pure classical gravity. D is a, differential 
operator. After a suitable number of integrations by part it is always possible 
to convert the quadratic part of the action to the form given in (1.2). We 
postpone discussion of possible boundary terms to the next sections. In the 
simplest case of quantum scalar field on the background of a classical geometry 
D is a, Laplacian with a mass term: 

D = Do -V^V^ + (1.4) 

Just this simple example is enough to illustrate the material of this section. 
Note, that in this case J has no contribution from the first variation of the 
classical action since gravity is not quantised. 
The path integral measure is defined by: 

1 = |p0exp(-(0,0)) . (1.5) 

Strictly speaking the right hand side of (1.5) is divergent. The essence of the 
condition (1.5) is that this divergence does not depend on external sources and 
on the background geometry and, therefore, may be absorbed in an irrelevant 
normalisation constant. The Gaussian integral (1.1) can be evaluated giving 

Z[J] = e-^='det-5(D)exp Qj^'V^ . (1.6) 

We stress that the equation (1.6) is true only if the operator D is self- 
adjoint. This means that D is symmetric or formally self-adjoint with respect 
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to the scalar product (.,.): 



((/.l,D(/.2) = Wl,02) (1.7) 

for any 0i, 02, and that the domains of definition of D and its adjoint coincide. 
We will not care about the second requirement since it involves mathemati- 
cal machinery [364] going beyond the scope of the present report. The first 
requirement (1.7) poses important restrictions on admissible boundary condi- 
tions. 

To become convinced that (1.7) is really necessary, one can calculate a 
finite-dimensional Gaussian integral with 




(1.8) 



and a, b, c and d being real constants first by completing the squares in 
the exponent, and then compare the result with det D (of course, one should 
remember a factor of tt) . The two results coincide if 6 = c. 

Let us return to the generating functional (1.6). To analyse the two mul- 
tiplies on the right hand side of (1.6) which depend on the operator D it is 
convenient to introduce the heat kernel 

K{t;x,y;D)^{x\eM-tD)\y). (1.9) 

This somewhat formal expression means that K{t; x, y; D) should satisfy the 
heat conduction equation 

{dt + D,)K{t-x,y-D) = Q (1.10) 

with the initial condition 

K{Q-x,y-D)=5{x,y). (1.11) 
For D — Do (1.4) on a fiat manifold M — M.^ the heat kernel reads: 

K{t;x,y;Do) = (47ri)-"/^ exp (^J^^^ - tm'^ . (1.12) 

The equations (1.10) and (1.11) can be checked straightforwardly. Let us con- 
sider a more general operator D which contains also a potential term or a 
gauge field. Then (1.12) still describes the leading singularity in the heat ker- 
nel as t — > 0. The subleading terms have a form of the power-law corrections ^ : 



^ On a curved space even the leading term must be modified (cf sec. 4.3). On 
manifolds with boundaries also half-integer powers of t appear in the expansion, 
and, consequently, 62^+1 7^ 0. 
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K{t; X, y- D) = K{t; x, y; Do) (l + thix, y) + eh^{x, y) + . . .) . (1.13) 

The coefficients hk{x, y) are regular in the hmit y ^ x. They are called the heat 
kernel coefficients^ . At coinciding arguments bk{x,x) are local polynomials of 
background fields and their derivatives. 

The propagator D^^{x,y) can be defined through the heat kernel by the 
integral representation 

roo 

D-\x,y)^ / dtK{t-x,y-D), (1.14) 

which follows from (1.9) if we suppose that the heat kernel vanishes sufficiently 
fast as t — >• cxD. We can formally integrate the expansion (1.13) to obtain 

D~\x,y) ^ 2(47r)-"/^ (^^^j i^-in+,+i(k " 2/|rn)^(a:, y) , 

(1.15) 

where ho = 1. By examining the behaviour of the Bessel function K^{z) for 
small argument z [251] we conclude that the singularities in the propagator at 
coinciding points are described by the first several heat kernel coefficients bk- 
Let us consider the part of the generating functional (1.6) which contains 
det{D). The functional 

W ^^\ndet{D) (1.16) 

is called the one-loop effective action. It describes the quantum effects due to 
the background fields in the one-loop approximation of quantum fields theory. 

To relate W to the heat kernel we shall use the arguments of [69]. 

For each positive eigenvalue A of the operator D we may write an identity 



InA = - 




(1.17) 



This identity is "correct" up to an infinite constant, which does not depend on 
A and, therefore, may be ignored in what follows'^ . Now we use lndet(£') = 
Trln(D) and extend (1.17) to the whole operator D to obtain 

1 /-oo A+ 

where 

K{t, D) = Tr (e-*^) = J d''x^K{t] x, x; D) . (1.19) 

^ We shall mostly use the coefficients which differ from by a normalisation 
factor. 

^ To "prove" this statement one has to differentiate both sides of eq. (1.17) with 
respect to A. 
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Here we have only presented some heuristic arguments in favour of eq. 
(1.18). A more rigorous treatment of functional determinants can be found in 
sec. 2.2. 

The integral in (1.18) may be divergent at both limits. Divergences at 
t — oo are caused by zero or negative eigenvalues of D. These are the infra 
red divergences. They will not be discussed in this section. We simply suppose 
that the mass m is sufficiently large to make the integral (1-18) convergent at 
the upper limit. Divergences at the lower limit cannot be removed in such a 
way. Let us introduce a cut off at t = A~^. 

It is now easy to calculate the part of W\ which diverges in the limit A — > oo: 
lydiv ^ _(47r)-V2y fin^^ 



2U+i)<n n-2j- 21 



+ H^){-myi\b2j{x,x) + 0{A°) 

2{j+l)=n 



(1.21) 



We see that the ultra violet divergences in the one-loop effective action are 
defined by the heat kernel coefficients bk{x,x) with k <n. 

On non-compact manifolds the integral of bQ{x,x) is divergent. This di- 
vergence is removed by subtracting a "reference" heat kernel (see sec. 6.1). 

Contributions from higher heat kernel coefficients bk, k > nio the effective 
action are not divergent and can be easily calculated yielding in the limit 
A ^ oo 

-hA^y^/^m- [ d-x^g E ^^^r(2j - n) . (1.22) 

^ 2j>n ^ 

This is nothing else than the large mass expansion of the effective action. This 
expansion is valid for relatively weak and slowly varying background fields. 
We have seen that the heat kernel expansion describes 

• short- distance behaviour of the propagator; 

• one-loop divergences and counterterms; 

• 1/m expansion of the effective action. 

We shall see below that heat kernel provides a natural framework for studying 

• quantum anomalies (sec. 7); 

• various perturbative expansions of the effective action (sec. 8); 

• selected non-perturbative relations for the effective action (sec. 9). 

Of course, in all these applications the heat kernel methods have to com- 
pete with other techniques. The main advantage of the heat kernel is that it 
delivers necessary information in terms of just few geometric invariants. This 
method does not make distinctions between different spins, gauge groups, etc. 
Even dependence of the space-time dimensionality is in most cases trivial. 
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Therefore, on one hand, just a single calculation serves then in many applica- 
tions. On the other hand, calculations in simple particular cases give valuable 
information on the general structure of the heat kernel. This property is espe- 
cially useful when one deals with complicated geometries (like in the presence 
of boundaries or singularities). During the last decade, many models which 
lead to such complicated geometries were very actively studied in theoreti- 
cal physics. The Dirichlct branes and the brane world scenario are the most 
popular but not the only examples. 

We have to mention also the limitations of the heat kernel formalism. It 
works less effectively in the presence of spinorial background fields, i.e. when 
there is mixing between bosonic and fermionic quantum fields. This problem 
is probably of the technical nature, so that the corresponding formalism may 
be developed some time in the future. A more serious drawback is that the 
heat kernel expansion is not applicable beyond the one-loop approximation. It 
is not clear whether necessary generalisations to higher loop could be achieved 
at all. 

It is not possible to write a review paper on heat kernel which would 
be complete in all respects, especially in the bibliography. A more compre- 
hensive treatment of many mathematical problems related to the heat kernel 
expansion can be found in [245,297]. The book by Kirsten [297] considers also 
specific physical applications as the Casimir energy and the Bose-Einstein 
condensation. The recent review paper [79] is devoted to the Casimir effect 
(see also [322]). The monographs by Birrell and Davies [69] and Fulling [224] 
remain standard sources on quantum field theory in curved space. Quantiza- 
tion of gauge theories is explained in [200]. In [107, 198, 189, 31, 191] the heat 
kernel expansion is treated from the point of view of quantum gravity and 
quantum cosmology. Useful information about properties of the zeta function 
can be found in [178]. One may also consult [182]. The DcWitt approach to 
the heat kernel and its' generalisations are described in [55]. The path integral 
point of view on the heat kernel can be found in [367] . 

Our primary goal is local formulae for the heat kernel coefficients. There- 
fore, in some cases global aspects will be somewhat neglected. 

This report is organised as follows. The next section contains necessary 
preliminary information on spectral geometry and differential geometry. There 
we define main geometric characteristics of the manifold and of the bound- 
ary. We discuss the zeta function (which defines the effective action) and the 
resolvent (which is a generalisation of the propagator). We relate the asymp- 
totics of these functions to the heat kernel coefficients. Through this report 
we work on Euclidean manifolds. A short remark on the analytical continua- 
tion to the Lorentzian signature is given at the end of section 2. In section 3 
we consider the most widely used models of quantum field theory and open 
bosonic strings. The one-loop dynamics in each of these models is defined by 
a second order differential operator which depends on an effective connection 
and on a matrix- valued potential. The connection and the potential serve as a 
basis of an invariant description of all that models in the language of spectral 
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geometry. These quantities are written down explicitly for each model. We 
also define suitable boundary conditions. In section 4 we consider the heat 
kernel expansion on manifolds without boundary. We introduce a simple and 
very powerful method and illustrate it by calculating several leading terms in 
the heat kernel expansion. We also briefly discuss some other methods and 
non-minimal operators. Section 5 is devoted to manifolds with boundary. The 
heat kernel expansion for standard Dirichlct, Neumann, and mixed boundary 
conditions is considered in some detail. We also describe less known oblique 
and spectral boundary conditions (these are the ones which contain tangential 
derivatives on the boundary or non-local projectors). We discuss loss of the 
so-called strong ellipticity for oblique boundary conditions which corresponds 
to the critical value of the electric field in string physics. With the example of 
spectral boundary conditions we illustrate appearance of non-standard (Int) 
asymptotics in the heat kernel. The results of sections 4 and 5 are valid on 
smooth manifolds with smooth potentials and gauge field. In section 6 we 
consider the case when either the background fields or the manifold itself have 
singularities. In particular, conical and domain wall singularities are consid- 
ered. 

Sections 7-9 are devoted to applications. In section 7 we relate confor- 
mal and chiral anomalies to certain heat kernel coefficients and re-derive the 
anomalies in several particular models. This is a very spectacular but also 
rather well known application of the heat kernel expansion. In section 8 we 
go beyond the power series in the proper time t. We consider mainly two par- 
ticular cases. The first one is the case when derivatives are more important 
than the potentials (so that we sum up leading non- localities) . The second one 
is the so-called low energy expansion which neglects derivatives of the back- 
ground fields starting with certain order (but treats all powers of background 
potentials and curvatures exactly). In that section we also review some results 
on the heat kernel on homogeneous spaces where spectrum of relevant oper- 
ators may be found exactly. In section 9 we consider two examples when the 
heat kernel can be used to obtain exact results for the effective action. The 
first one is the famous Polyakov action. The second example is exact relations 
between the effective actions in dual theories. In this section (in contrast to 
section 8) we don't have to neglect derivatives or powers of the background 
fields. Section 10 contains concluding remarks. 
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2 Spectral functions: heat kernel, zeta function, resolvent 

2. 1 Differential geometry and the operators of Laplace type 

Let M be a smooth compact Riemannian manifold ^ of dimension n with 

smooth boundary dM. Wc shall also consider the case when the boundary 
dM is empty. Let V be a vector bundle^ over M. This means that there is a 
vector space attached to each point of the manifold. For example, this could 
be a representation space of a gauge group or of the space-time symmetry 
group. Sections of V are smooth functions bearing a discrete index which 
could correspond to internal or spin degrees of freedom. 

We study differential operators on V. We restrict ourselves to second order 
operators of the Laplace type. Locally such operators can be represented as: 

D^-ig'^'^d^d^ + a'^d. + b) (2.1) 

where g'^'^ is the inverse metric tensor on M; a'^ and b are matrix valued func- 
tions on M. There is a unique connection on V and a unique endomorphism 
E oiV (another matrix valued function) so that 

D = -{g^^V^V. + E), (2.2) 

where the covariant derivative V = V^^^ +u contains both Riemann V'^' and 
"gauge" (bundle) cu parts. We may express: 

^5 = ^^?.5(a" + ^''"r^/Jy) and (2.3) 
E = b- g'^'id^uj, + uj^uj^ - a;,r,/) , (2.4) 

where 

r./ = g'"'\ {d^^gup + d^g^, - d,g^,) (2.5) 

is the Christoffel symbol. ly is the unit operator on V. 

It is important to understand how general the notion of the Laplace type 
operator is. The most obvious restriction on D is that it is of second or- 
der, i.e. it contains second derivatives, but does not contain higher derivative 

^ This simply means that we assume that there is a positive definite metric tensor 
gij,^ on M. 

^ Many physicists strongly dislike vector bundles. Nevertheless, there are two good 
reasons for using the fibre bundles in this paper (in parallel with more familiar nota- 
tions of matrix- valued functions, gauge fields, etc). First, our simplifying comments 
and examples may help the reader to understand mathematical literature on the 
subject. Second, one of the main ideas of this report is to reveal some universal 
structures behind the heat kernel expansion. In particular, we shall see that there 
is not much difference between different spins and symmetry groups. The vector 
bundle language seems to be the most adequate language to achieve this goal. The 
reader may consult the excellent review paper by Eguchi, Gilkey and Hanson [177]. 
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parts. In this paper we shall also consider first order operators (Dirac oper- 
ators, for example). Second, the operator D is a partial differential operator. 
This excludes negative or fractional powers of the derivatives. The operators 
containing such structures are called pseudo-differential operators. More infor- 
mation on spectral theory for these operators can be found in [383,258,259]. 
Third, the operator (2.1) has a scalar principal part. This means that the sec- 
ond derivatives in (2.1) are contracted with the metric, and the internal index 
structure of the second derivative term is trivial. Such operators are also called 
minimal. Non-minimal operators will be briefly considered in sec 4.4. 
We can define local invariants associated with cu and g. Let 

= a.r;;^ - a.r;;, + r^/^, - r^^r^^ (2.6) 

be the Riemann curvature tensor, let i?^^ := ^va be the Ricci tensor, and let 
R := be the scalar curvature. With our sign conventions, 7? = 2 on the unit 
sphere S'^ in the Euclidean space. Let Roman indices i, j. A;, and / range from 1 
through the dimension n of the manifold and index a local orthonormal frame 
(vielbein) {ei,...,e„} for the tangent space of the manifold. In components 
we have: e^e^g'^,^ = 5jk, e^e^5^^ — g'^^ . The inverse vielbein is defined by the 
relation e^^e^ = 5^. These two objects, and e^, will be used to transform 
"curved" indices (;U, z/, p) to "flat" ones (i, j, k) and back. In Euclidean space 
there is no distinction between upper and lower flat indices. 

As usual, the Riemann part of the covariant derivatives contains the 
Christoffel connection so that 

Vjfl^. = d,v^ - V^^v, (2.7) 

for an arbitrary vector v^,. To extend this derivative to the objects with flat 
indices one has to use the spin-connection o"^: 

V^v^ = d^v^ +<jl^Vk. (2.8) 

The condition Vj^e^ = yields 

< = err^^ej - e^d.el . (2.9) 

Let 0,^1, be the held strength of the connection uj: 

fl^iu = d^uj„ - d^uj^ + uj^uj^ - uj^uj^. (2.10) 

The covariant derivative V acts on both space-time and internal indices. For 
example, 

v,fi^, = d,n^, - T;^n^, - v^^n^^ + [u^, n^,] . (2.11) 

If the boundary DM is non-empty we have more invariants. Let e„ be 
inward pointing unit vector fleld. Let Roman indices a, 6, c and d range from 
1 to n — 1 and index a local orthonormal frame for the tangent bundle of dM. 
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Let Lab '■— ^ab ^he second fundamental form (extrinsic curvature) of the 
boundary. We use the Levi-Civita (spin) connections and the connection u 
to covariantly differentiate tensors of all types. Let ';' denote multiple covari- 
ant differentiation with respect to the Levi-Civita connection of M and let 
':' denote multiple tangential covariant differentiation on the boundary with 
respect to the Levi-Civita connection of the boundary; the difference between 
';' and ':' is measured by the second fundamental form. Thus, for example, 
E-a = E-a since there are no tangential indices in E to be differentiated. On 
the other hand, E-ab 7^ E,ab since the index a is also being differentiated. More 
precisely E.^b — E-.ab ~ LabE-^. Since L is only defined on the boundary, this 
tensor can only be differentiated tangentially. 

Consider an example of the circle in the plane M^. The line element 
has the form 

{dsf = {drf + r\de)\ (2.12) 

where 0<r<cxD, 0<6'< 27t. Then gi^^ = diag(l,r^). is defined by the 
condition r = tq. We may choose — r, — —1 (the minus sign appears 
because is an inward pointing unit vector for the disk with the boundary 
S^). Then the second fundamental form of is 

Lu^eyA^e (2.13) 

r=ro Tq 

In general, on S^~^ considered as a boundary of the ball in M** the extrinsic 
curvature is Lab = ^Sab- 

If the boundary dM is non-empty, one has to define boundary conditions 
for the field 0. A convenient way to write them down is: 

B(j)^0 (2.14) 

where B is called the boundary operator. In general, the operator B calculates 
a linear combination of the boundary data for any given function 0. If D is of 
Laplace type, the boundary data include value of the function at the boundary 
and value of it's first normal derivative. The most frequently used choices are 
the Dirichlet and Neumann boundary operators which we denote B^ and B~^ 
respectively: 

i3-0 = 0|aM, (2.15) 
-B+0=(0;„ + 50)|9M, (2.16) 

where S is a matrix valued function defined on dM. The boundary conditions 
(2.16) are also called Robin or generahsed Neumann. In some literature they 

are called mixed boundary conditions. We shall not use this latter terminology. 
The name "mixed" is reserved for another type of boundary conditions. Let 
n~ and n+ be two complementary projectors defined on V\qm, O^^^Y = n^, 
n+ + n- = I. There is a decomposition T^Iqm — Vn ® Vd, where Vn d — 
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n y laM- Decompose also = ^ © and set 

i30:=0^©(0;^ + 50^)|9M. (2.17) 

The matrix valued function (endomorphism) 5" acts on Vn only, 5" = U'^S — 
SH'^. In other words, we define Dirichlet boundary conditions on Vd and gen- 
eralised Neumann boundary conditions on Vn- For obvious reason the bound- 
ary conditions (2.17) will be called mixed. In sec. 3 we shall see that natural 
boundary conditions for spinor and vector fields are of this type. 

2.2 Spectral functions 

For the boundary conditions considered in this section as well as on man- 
ifolds without a boundary the operator exp{—tD) with positive t is trace class 
on the space of square integrable functions L'^{V). This means that for an 
auxiliary smooth function f on M 

K{t,f,D) = TrL2{fexp{-tD)) (2.18) 

is well defined. We also write 

K{t, f,D)^ [ (Tx^iryKit- X, x; D)f{x) , (2.19) 

JM 

where K{t] x, x; D) is any x limit of the fundamental solution K{t; x, y; D) 
of the heat equation (1.10) with the initial condition (1.11). If there is a 
boundary, the kernel K{t] x, y; D) should also satisfy some boundary condi- 
tions BxK{t; X, y;D) = in one of the arguments. We stress that K{t; x, y; D) 
is a matrix in the internal indices, try denotes the trace over these indices. 

Let D be self-adjoint. This implies that in a suitable basis in the internal 
space the matrix o;^ is anti-hermitian and E is hermitian. Let {4>\} be a 
complete basis of orthonormal eigenfunctions of the operator D corresponding 
to the eigenvalues {A}. Then 

Kit; X, y;D) = J2 0i(x)0A(?/)e-*' . (2.20) 

A 

We shall almost exclusively work either on manifolds without boundaries, 
or on manifolds with boundaries with the fields subject to local boundary 
conditions (2.15), (2.16), or (2.17). In all these cases there is an asymptotic 
expansion as i i [381, 382, 253] ^ : 

TrL2{fexp{-tD)) = ^ t(^-")/V(/, ^) • (2.21) 

jfc>0 

^ These papers contain also a method allowing to calculate the coefficients of the 
expansion. The method is, however, too complicated to use it on practice. 
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This expansion is valid for almost all boundary conditions appearing in appli- 
cations to physics. There are, however, some exceptions which will be discussed 
in sec. 5.4. 

The coefficients and the coefficients bk introduced in the previous sec- 
tion (cf. (1.13)) are related by the equation (m — for the simplicity): 

Ukif, D) = (47r)-"/2 / (Tx^hkix, x)f{x) . (2.22) 

Note, that the definition (1.13) is valid on fiat manifolds without boundary 
only (although generalisations to other cases are possible) . 

The key property of the heat kernel coefficients ak is that they are locally 

computable in most of the cases. This means that they can be expressed in 
term of the volume and boundary integrals of local invariants. 

For a positive operator D one can define the zeta function [325, 383] by 
the equation: 

C(s,/,L>) = IV^2(/L>-^). (2.23) 
The zeta function is related to the heat kernel by the integral transformation 

oo 

as,f,D)^r(s)-' J dtf-'K(t,f,D). (2.24) 



This relation can be inverted, 

K{t, LD)^^Jds rT(s)C(s, /, D) , (2.25) 

where the integration contour encircles all poles of the integrand. Residues at 
the poles can be related to the heat kernel coefficients: 

akif, D) = Res,=(„_fc)/2 (r(s)C(s, f,D)) . (2.26) 

In particular, 

a„(/,D)=C(0,/,I^). (2.27) 
Z eta-functions can be used to regularize the effective action [165,274]. 
The regularization is achieved by shifting the power of t in (1.18): 

^K{t,D), (2.28) 

where fx is a constant of the dimension of mass introduced to keep proper 
dimension of the effective action. The regularization is removed in the limit 
s — > 0. Eq. (2.28) can be considered as a definition of the regularized effective 
action without any reference to (1.18). One can also rewrite (2.28) in terms of 
the zeta function: 



Ws^-lfi'ns)C{s,D), (2.29) 



where C{s,D):=C{s,l,D) 
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The gamma function has a simple pole at s = 0: 

r(s) = --7E + C»(s), (2.30) 

S 

where je is the Euler constant. The regularized effective action (2.29) has also 
a pole at s = 0: 

= - 7s + In/i') C(0, D) - ic'(0, D) . (2.31) 

According to (2.27) the divergent term in the zcta function regularization is 
proportional to an{D) (cf. (1-21) for another regularization scheme). The pole 
term in (2.31) has to be removed by the renormalization. The remaining part 
of at s = is the renormahsed effective action: 

ly- ^ _ 1 ^.(0^ ^) _ 1 in(^2^C(0, D) , (2.32) 

where we have introduced a rescaled parameter fi^ = e^^'^ji^. In this approach 
n"^ describes the renormalization ambiguity which must be fixed by a suitable 
normalisation condition. Let us remind that here we are working on a com- 
pact manifold. On non-compact manifolds ({s, D) may have divergent contri- 
butions proportional to the volume. Such divergences are usually removed by 
the subtraction of a "reference" heat kernel (see sec. 6.1). 

Together with (1.16) equation (2.32) yields a definition of the functional 
determinant for a positive elliptic second order operator which is frequently 
used in mathematics: 

Indet(D) = -C'(0,D) -ln(/.2)C(0,D) . (2.33) 

Note, that the definitions (2.23), (2.33) are vahd for positive operators 
only. Elliptic 2nd order differential operators have at most finite number of 
zero and negative modes ^ which must be treated separately. However, one can 
extend the definition of the zeta function to operators with negative modes: 

C(s,L>) = $:|A|-% (2.34) 

where the sum extends over all non-zero eigenvalues A. One can also define 
another spectral function in a similar way: 

r){s,D)^Y.^\gn{X)\Xr. (2.35) 

This function is especially useful in spectral theory of Dirac type operators 
where r](0,p) measures asymmetry of the spectrum. 

^ On manifolds with boundaries one has to require strong ellipticity of the boundary 
value problem (cf. sec. 5.4). 
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Another function which is frequently used especially in the mathematical 
literature is the resolvent (or, more precisely, its' powers): 



R'(z) := (D + z')-'. (2.36) 

If D is on operator of Laplace type subject to "good" boundary conditions^ 
and if / sufficiently large (/ > n/2) there is a full asymptotic expansion 

Tr(R'(^)) = ^ U ,{D)z-'^+--' (2.37) 

as ^ — > cxD. The coefficients are the same as in the heat kernel expansion 
(2.21). 

2.3 Lorentzian signature 

Locality of the heat kernel coefficients in the Euclidean domain can be 
easily understood by examining the free heat kernel (1.12). For small t the 
first term in the exponential strongly suppresses non-local contributions. For 
Lorentzian metrics the squared distance function {x — yY is no longer positive 
definite. Therefore, the simple arguments given above do not work. A partial 
solution to this problem is to consider the "Schrodinger" equation 

{idr- D)K{T;x,y;D)^0 (2.38) 

for the kernel K instead of the heat conduction equation (1.10). Then K 
oscillates at large distances. However, even though non-local contributions 
to K oscillate furiously as r — > they are not small. Consequently, local 
asymptotic series do not exist in many cases. A discussion on this point can 
be found in chapter 9 of [224]. 

Of course, the heat kernel expansion can be used also on Lorentzian man- 
ifolds at least for the renormalization theory where non-local terms are of less 
importance. Counterterms are still defined by the same heat kernel coefficients 
with the same functional dependence on the Lorentzian metric. Explicit def- 
initions with the imaginary "proper time" r can be found in [147,69]. One 
should note that some background fields receive an imaginary phase when 
being continued to the Euclidean domain (cf. sees. 3.2 and 3.3). 



A more precise meaning of this restriction will be explained in sec. 5.4. 
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3 Relevant operators and boundeiry conditions 



The operator of Laplace type is not necessarily the scalar Laplacian. In 
fact, in almost all models of quantum field theory the one loop effective action 
is defined by an operator of this type. This can be demonstrated by bringing 
relevant operators to the canonical form (2.2). In this section wc give ex- 
plicit construction of the connection uo and the endomorphism (matrix-valued 
potential) E for scalar, spinor, vector and graviton fields. We also describe 
appropriate boundary conditions. 



3.1 Scalar fields 



Consider first the example of the multi-component real scalar field in 
n dimensions. The action reads 

£ = / (rx^{g^''V^^^VM + U{^)^im^^^), (3.1) 

J M 

where ^ is the conformal coupling parameter, U is a. potential. The covariant 
derivative contains the background gauge field G^^: V^^^"^ = d^^"^ -\- 
G^^^^ . Gfj_ is antisymmetric in internal indices A, B. To evaluate the one- 
loop efi'ective action one should expand (3.1) around a background field 
$ = $ + 0, and keep the terms quadratic in fiuctuations: 

A = (T'x^cp^ (-(V^V^)-^^ + \{U{^)'Y'' + iR5'''') cp"" 

- [ (P-^xVh(j)^Vn(p^, (3.2) 

JdM 

h is the determinant of the induced metric on the boundary. The inner product 
for quantum fields reads 

<</)!, 02 >= / (rx^(l>^(t>^ . (3.3) 

JM 

The operator D is defined by the bulk part of the action (3.2): 

-(V'^V^)'^^ + -(^7($)'y^ + iR5^^ . (3.4) 
For a special choice of the parameter ^ : 

the operator D (3.4) is conformally covariant (if also U{^)" = 0). To bring 
the operator (3.4) to the canonical form (2.2) we introduce 

^AB ^ ^AB ^ ^AB ^ _l(C/($)")^5 _ ^R5^B (3 g) 
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For this case, Q.^^, is just the ordinary Yang-Mills field strength. 

The operator (3.4) is symmetric with respect to the inner product (3.3) 
if the surface integral 

/ r-'x Vh{4>^Vn(t>^ - ct>^VnCl>t) (3.7) 
JdM 

vanishes for arbitrary 0i and 02 belonging to its' domain of definition. This 
may be achieved if one imposes either Dirichlet 

(t>\M = (3.8) 

or modified Neumann 

(V„0^ + 5^^0^)|aM = O (3.9) 

boundary conditions. S"^^ is an arbitrary symmetric matrix. Note, that the 
integral (3.7) vanishes also if S^^ is an arbitrary symmetric (differential) op- 
erator on the boundary. 

For the Dirichlet conditions (3.8) the boundary term in (3.2) vanishes 
automatically. To ensure absence of the surface term for the modified Neumann 
conditions (3.9) one should add to (3.1) an appropriate surface action. 

3.2 Bosonic string 

Our next example is the non-linear sigma model in two dimensions de- 
scribed by the action 

£M = d^x (ygGAB{X)g^'d^X^d,X'' + e^'' BAB{X)d^X^d,X'') 

+ / AedX^ . (3.10) 

JdM 

From the point of view of two-dimensional world sheet the fields X'^i^x) are 
scalars. In string theory they arc interpreted as coordinates on a rf-dimensional 
target manifold with the metric Gab{X). e^^ is the Levi-Civita tensor density, 
e^^ — — e^^ = 1. Bab{X) is an antisymmetric tensor field on the target space. 
Ab{X) is the electromagnetic vector potential. The action (3.10) describes 
charged open stings. For simplicity, we absorb the inverse string tension a' 
into a field redefinition. We do not include tachyon and dilaton couplings in 
the bulk or on the boundary. 

Usually the term with the S-field gets an imaginary coefficient in the 
Euclidean space. Since the physical space-time has Minkowski signature it is 
not especially significant which way of continuation to the Euclidean space has 
been chosen provided the results are properly continued back to Minkowski 
space after the calculations. As we will see below, real coefficient in front of 
the S-field leads to a well-defined spectral problem. This situation is in close 
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analogy with the continuation rules for the axial vector field in the spinor 
determinant [14, 13]. The other way to deal with the field B is to keep the 
coefficient of the i^-term in (3.10) imaginary at the expense of introducing a 
more sophisticated conjugation operation [352] containing the sign-reversion 
of the S-field. The same refers to the electromagnetic potential Ab- 

The field X enters the action (3.10) at many places making the back- 
ground field expansion a quite cumbersome procedure. The most economic 
way to arrange such an expansion and to calculate higher derivatives of the 
action (3.10) is to introduce the geodesic coordinates in the target space. A 
detailed explanation of the method as well as further references can be found 
in [87]. Consider the target space geodesies defined independently at each point 
of the two-dimensional world surface and parametrised by the arc length s in 
the target space. They satisfy the usual geodesic equation 

where ^bc^-^) the Christoffel connection corresponding to the target space 
metric Gab- Let us supplement the equation (3.11) by the initial conditions 

X\x, 0) = X\x) , j-X\x, 0) = , (3.12) 

where X is the background field. parametrises deviations from X and, 
therefore, can be identified with quantum fluctuations. The fc-th order term 
of the expansion of the action (3.10) around the background X is given by 

^^ = ^^A^(^))l-o. (3.13) 

Higher order derivatives oi X{x, s) with respect to s can be traded for the first 
derivatives by means of the geodesic equation (3.11) and then replaced by ^ 
at s = with the help of the initial conditions (3.12). 

It is convenient to introduce the Riemann curvature of the target space 
metric Gabi T^abcd, and a 3-index field strength 

Habc = QaBbc + dsBcA + dcBsA ■ (3-14) 
The covariant derivative V is 

"^.e = + lic{X){d,X<')e + \e,^{d,X'')H^Bcf • (3.15) 
The quadratic part of the action (3.10) reads: 
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M 



~{d,X<'){d,X'')e^^eeDAHBCD 

-r\{d,X^){d^X^)HABcH'' DE^e) (3.16) 
+ / iD,ee{BAB - Fab) + {d,X^)DA{BBc - Fbc)^^^) , 

JdM ^ ' 

where r is the arc length along the boundary. Fbc = dsAc — OcAb- The 
covariant derivatives D^- and Da contain the Christoffel connection on the 
target space (but not Habc as the full covariant derivative V, eq. (3.15)). 
The natural inner product in the space of fluctuations ^ reads 

< C(i),C(2) >= j^d''x^GAB{X{x))i^,){x)if^){x) . (3.17) 

The volume part of the action (3.16) has now the canonical form D^) with 

the operator D (2.2) which is obviously of Laplace type. The connection uj is 
defined in (3.15) and the endomorphism E can be easily extracted from (3.16). 
The operator D is formally self-adjoint with respect to (3.17) if we impose the 
boundary conditions of Neumann type 

B^^{Vne + S^^^)\9M^0 (3.18) 

with arbitrary operator S which should be symmetric with respect to the 
restriction of (3.17) to the boundary: 

cir(e(t)5|e(2)s - ^^2)S^^ii)B) = (3.19) 

There is a preferable choice of the boundary operator. Let us vary the action 
(3.16) with respect to the fluctuation field ^: 

ldS2= [ d'x^/h{dOD^- [ dT{5i)Bi. (3.20) 

Z JM JdM 

Now we require that the boundary integral in (3.20) vanishes for arbitrary 5^. 
Hence we arrive at the boundary conditions (3.18) with the operator S given 

by 

= ^(rv. + v.r)^ + 5^ , = - Fa^ , (3.21) 

-S"^ = -^{dfiX'^) H^^c{Bdb — Fdb) + bc{Bd^ — F^^) 

+^(a,X^) [l>^(Bbc - Fbc) + Db{B^c - F^c) 

Note that the operator in (3.21) is not of the ordinary Neumann (or Robin) 
type since it contains tangential derivatives on the boundary (cf. sec. 5.4.1). 
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The variation (3.20) vanishes also if we choose Dirichlet boundary con- 
ditions for some of the coordinates of the string endpoints. Namely, we can 
take a projector n_|_ and impose (3.18) on (n_|_)^^'' and {1 — II+)^S^^\qm = 0. 
Physically this means that the endpoints of the bosonic string are confined 
to a submanifold in the target space. Such configurations [134] are called the 
Dirichlet branes. 

3.3 Spinor fields 

The action for the spinor fields ■0 

£ - / d^'xy/gippilj (3.22) 

contains a first order operator p of Dirac type. In Euclidean space the conju- 
gate spinor 'ijj is just the hermitian conjugate of ■0: t/J = By definition, an 
operator p is of Dirac type if its square D — p^ is of Laplace type. Spectral 
theory of general operators of Dirac type both on manifolds without bound- 
aries and with local boundary conditions on manifolds with boundaries can 
be found in [88,89]. Here we consider some physically motivated particular 
cases only. Let us introduce the Euclidean Dirac 7-matrices which satisfy the 
Clifford commutation relations: 

lulu + lulv = 2gf,„ . (3.23) 

The 7-matrices defined in this way are hermitian, ■jj^ = 7^. We also need the 
chirality matrix which will be denoted 7^ independently of the dimension. It 
satisfies 

= 7^7,. = -7m7' • (3.24) 

Prom now on we suppose that the dimension n is even. We fix the sign in 7^ 
by choosing 

_ n(n — 1) 

7' = ^-e'^-'^7M---7.- (3.25) 

Let p be the standard Dirac operator in curved space with gauge and axial 
gauge connection 

P = {d, + ^ [7., 7pKf + A, + z47^) . (3.26) 

Here a'j^'' is the spin- connection (2.9). and are vector and axial vector 
fields respectively taken in some representation of the gauge group. Both A^ 
and are antihermitian in the gauge indices. The operator p (3.26) is for- 
mally self-adjoint in the bulk. The operator D = p^ is of Laplace type (2.2) 
with 10 

The present author is grateful to Valery Marachevsky for his help in deriving and 
checking eqs. (3.27), (3.28) and (3.37), see also [136]. 
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'^M = ^[7., lM7 + + ^[7^. 7.]^'"7' , 

-^(n-3)[7^7l[4'^'] (3-27) 

with obvious notations F^^ = d^A^-d„A^ + [Af„ A^], D^Al = d^Al-VP^^Al + 
[An,A^J[. The expression for Q,^^, is a httle bit lengthy: 

= F^. - [A% Al] - lrYRa,,u - ii'Yil.D^Al - ^,D,Al) 
Mi'A% + [4, 417^7, - [Al Al]Yl, 

-YAl^.^Al^, + 7M^7,7M^7^ , (3.28) 

where 

Al = d,Al - d^Al + [A„ Al] - [A,, Al] (3.29) 

Consider now manifolds with boundary. The specific feature of the action 
(3.22) is that it contains first order derivatives only. Consequently, boundary 
conditions should be imposed on a half of the spinor components. Let these 
be Dirichlet boundary conditions 

n.V^IsM = (3.30) 
where n_ is a hermitian projector, = n_, lit = n_. Due to the hermiticity 

V'n_|aM = (3.31) 
Following Luckock [308] let us consider a family of the projectors 

n_ = ^(l + 7-exp(zg75)) , (3.32) 

where is a scalar which can depend on the coordinate on the boundary. 

To make the operator p formally self-adjoint including the boundary we 
must require that 

/ r-^x v^V'i7nV'2 = (3.33) 

JdM 

for all -01^2 satisfying the boundary conditions (3.30). Since ■^i and -02 are 
arbitrary, the projector n_ should satisfy 

(i-n_)7"(i-n_) = 0. (3.34) 

This condition yields q — ±7r/2. The projector (3.32) takes the form 

n_ = ^ (1 ± i^^Y) . (3.35) 



23 



To formulate the spectral problem for the second order operator D — 
we need boundary conditions for the second half of the spinor components. 
The relevant functional space should be spanned by the eigenfunctions of the 
Dirac operator^. It is clear that on this space the functions ^t/) should satisfy 
the same boundary conditions (3.30) as the ■^'s themselves: 

njp^l:\9M^Q. (3.36) 

Let us adopt the choice (3.35) for n_. By commuting n_ with^ in (3.36) we 
obtain 

(v„ + 5)n+v^|sM = o, 5 = -^L„„n+, n+ = i-n_. (3.37) 

We remind that L^b is extrinsic curvature of the boundary. The boundary 
conditions (3.30), (3.37) with (3.35) are mixed (cf. (2.17)). Spectral geometry 
of the Dirac operator with these boundary conditions has been thoroughly 
studied by Branson and Gilkey [89]. 

We can generalise the boundary conditions presented above by consider- 
ing non-hermitian projectors, III ^ n_. Then the boundary condition for the 
conjugated spinors reads: 

V^nLlaM^O. (3.38) 
Instead of (3.34) we have the condition 

(i-nl)7"(i-n_) =0, (3.39) 

which yields 

n_ = ^ (l + i7"75e''(=^")^') (3.40) 

with an arbitrary real function (or even with an arbitrary hcrmitian matrix 
valued function) r{x'^). 

In the Minkowski signature space the spinor conjugation includes 7°. 
Therefore, the boundary conditions on the conjugate spinor (3.31) are changed. 
Roughly speaking, to continue n_ to the Minkowski signature space time one 
has to replace 7^7^ by 7" in (3.35) and (3.40) and to take into account powers 
of i which appear in the Dirac gamma matrices. 

In particle physics, interest to the boundary conditions defined by (3.35) 
is due to the MIT bag model of hadrons proposed in [124,125,143]. This model 
was modified later [126,111,95,397,365,250] to include a chiral phase on the 
boundary in a manner similar to (3.40) (for a review see [273]). Renormaliza- 
tion of quantum field theory with the boundary conditions defined by (3.35) 
was considered by Symanzik [394]. 

As we have already mentioned above, bag boundary conditions defined 
by the projector (3.35) are a particular case of mixed boundary conditions 
(2.17). They will be considered in detail in sec. 5.3 (see also sec. 5.1 for further 
references to calculations in a ball). Chiral bag boundary conditions defined 
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by (3.40) with r 7^ are considerably more complicated because the equa- 
tion (3.36) contains a mixture of normal and tangential derivatives. We refer 
to [281, 203, 204, 138, 421, 139, 137, 62, 199] for more calculations of spectral 
functions in this latter case. 

The boundary conditions considered in this section are local, i.e. they 
treat the fields at each point of the boundary independently. One can also 
define global houndaij conditions for the Dirac operator (cf. sec. 5.4.2). 

Due to the fermionic nature of the spinor field the path integral over tjj 
and ■0 gives determinant of the operator ^ to a positive power, 

J Vi^V'il^ex^i-Cii),^))) = det|> (3.41) 

where the action C is given by (3.22) and -0, -0 are complex Dirac spinors. 



3.4 Vector fields 



Consider the Yang-Mills action for the gauge field (Greek letters from 
the beginning of the alphabet label generators of the gauge group) : 

C = \jj^x^F^,F^''\ (3.42) 

where, as usual, F"^ = 9//^" — d^A"^ + c'^^A^^A^. c^^ are totally antisymmetric 
structure constants of the gauge group. Let us introduce the background field 
B^l by the shift A^ + A'^. From now on A^ plays the role of quantum 

fluctuations. The quadratic part of the action reads: 

C2 = - I (Tx^ [-a^'^v'^v^a;^ + a^v^vm;^ + a'^a^r^''' 

2 Jm l 

+2F;,(i?)c^^A^M^^; 

+J / (^-^x^^h[A'iV''A'^-A''''VnA'i\. (3.43) 

2 JdM 

The covariant derivative V contains both metric and gauge parts: V^A'^ — 
dyA'^ + B^c'^^A^ — TP^A^. One should impose a gauge condition on the fluc- 
tuation A^. We choose 

V'A^ = . (3.44) 

In the gauge (3.44) the bulk term in (3.43) defines an elliptic operator of the 
Laplace type with 

M7p = BX^K - , (3.45) 

(ET/, = -R% + '^nSVA. . (3.46) 



The field strength corresponding of the connection (3.45) reads 

(fi,.);^- = + F(S)J,c?/; . (3.47) 
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Note, that all the quantities above (3.45) - (3.47) are matrix functions with 
both gauge and vector indices. 

On a manifold without boundary the operator D with (3.45) and (3.46) 
is symmetric with respect to the standard inner product in the space of the 
vector fields 

< A(^) , >= y/gdTx ^(i)«^^A(^2)a _ (3_4g) 

The ghost operator corresponding to the gauge (3.44) is just the ordinary 
scalar Laplacian, D^^ — —'SJ\■9H^^SJ^9^\ with the connection 

The one-loop path integral reads: 

Z{B) = det(L>)7 det(L>[^'^])5 (3.50) 

where the first determinant is restricted to the fields satisfying the gauge 
condition (3.44). Note, that pure gauge fields A = are zero modes of the 
total operator D + VV of the bulk action (3.43) only on shell, i.e. when the 
background field B satisfies the classical equation of motion \/^F{B)'^^ = 0. 
Therefore, the path integral is gaugc-indcpcndcnt on shell, but it depends on 
the gauge choice off shell. For example, the Feynman gauge path integral 

ZpiB) = det(L')-^ det(i:>[»'*l) , (3.51) 

where the first determinant is calculated on the space of a// vector field, is equal 
to the Z{B) defined in (3.50) only on shell. However, physical predictions of 
the two path integrals are, of course, equivalent. 

The path integral (3.51) can be obtained by adding the gauge fixing term 

£ J_/ rx^{V,A^)' (3.52) 

zk. Jm 

to (3.43) with K — 1. The case k ^ 1 yields a non- minimal operator on the 
gauge field fluctuations (cf sec. 4.4). 

One can define "total" heat kernel coefficients for the path integral in a 
certain gauge x(^) 

= akiD^) - 2afc(Df 1) , (3.53) 

where is defined by the action (3.43) with the gauge fixing term x^iA), and 
the ghost operator is D^^^ = x(V{). Even on shell only the coefficient a^°* is 
gauge-independent. Only this coefficient contains information on the one- loop 
divergences in a gauge-invariant regularization like the C-function one. 

On a manifold with boundary one should impose boundary conditions 
on gauge fields and ghosts. The boundary conditions should be gauge invari- 
ant. Consider a more general set-up when we have some quantum fields $ 
and a linearised gauge transformation $ — > $ -|- 5^^ with local parameter ^. 
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Boundary operator B defines gauge invariant boundary conditions 



i3$ = (3.54) 
if there exist boundary conditions for the gauge parameter ^ 

= (3.55) 

such that 

B6!:(!> = . (3.56) 

The equation (3.56) means that the functional space defined by the boundary 
conditions (3.54) is invariant under the gauge transformations provided the 
gauge parameter ^ satisfies (3.55). Upon quantisation (3.55) become boundary 
conditions for the ghost fields. The condition (3.56) ensures validity of the 
Faddeev- Popov trick on a manifold with boundary and guarantees gauge- 
independence of the on-shell path integral [405]. 

Since the boundary term in (3.43) is diagonal in the gauge index a we can 
consider the case of an abelian gauge group and drop a from the notations. 
The general non-abelian case does not offer considerable comphcations. There 
are two sets of local boundary conditions which satisfy the gauge- invar iancc 
requirements above. The first set is called absolute boundary conditions. It 
reads 

■^n\dM — , dnAa\dM — (^n^ab — Lab)Ab\QM — , 

V„e|aM = 0. (3.57) 
The second set 

(V„ — Laa)An\dM — , ^al^M = , 

C\dM = (3.58) 

is called relative boundary conditions. The projectors on the Dirichlet and 
Neumann subspaces and the endomorphism S are = ^in^jm (n"'")^ = 

5ij ~ ^in^jn-i Sij = Lab^ia^jb and (IT )jj = 5ij SinSjn, (Tl^)ij = SinSjn, Sij 

I.. ' for the absolute and relative boundary conditions respectively. 
It is straightforward to check gauge invariance of the absolute boundary 
conditions (3.57). A bit more job is required to show that the condition (3.58) 
for the normal component is also gauge invariant. The key observation is 

that near the boundary the scalar Laplacian can be represented as DIs'i] _ 
~(V„— L(ia)V„+£^, where £ docs not contain normal derivatives and, therefore, 
£i\dM = under the boundary conditions (3.58) for the ^. Consequently, 

(V„ - Laa)S^An\aM = (V„ - La^)Vn^\9M = {-D^'^k + ^OldM = (3.59) 

where D^sh]^^ 

QM = on the eigenfunctions of the operator which is 

enough for our purposes. 
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The boundary term in the action (3.43) can be rewritten as 

I [ r-'x^[A,{A^.,, + L,,A, - A.„)] . (3.60) 

Z JdM 

This boundary action vanishes for both types of the boundary conditions 
(3.57) and (3.58). Hence the operator D is symmetric for these boundary 
conditions. Another remarkable fact is that the fields A-^ satisfying the gauge 
condition (3.44) are orthogonal to the gauge transformations. Indeed, 

< A^, V{ >= - / r-^xVhA^i = (3.61) 

JdM 

for both (3.57) and (3.58). 

We should also check whether the gauge condition (3.44) is indeed com- 
patible with the boundary conditions (3.58) and (3.57). For an arbitrary A 
there should exist unique ^ such that 

V^(A^ + V^O = 0. (3.62) 

We rewrite (3.62) as 

VM^ = -Ae . (3.63) 

Let us start with relative boundary conditions (3.58). In this case the left hand 
side of (3.63) satisfies Dirichlet boundary conditions. The scalar Laplacian A 
for these boundary conditions is invertible. Therefore, a solution for (3.63) 
always exists, and it is unique. The case of absolute boundary conditions is a 
bit more involved. One should take care of a zero mode in the ghost sector. 
We leave this case as an exercise. A more extensive discussion of compatibility 
of gauge and boundary conditions can be found in [405]. 

3. 5 Graviton 

We start this section with the Einstein-Hilbert action on a four dimen- 
sional Euclidean manifold without boundary: 

where R is the scalar curvature. Gat is the Newton constant, A is the cosmolog- 
ical constant. As usual, let us shift the metric gf^^ — > g^v+hni,. From now on g^i, 
will denote the background metric, h^^ will be the quantum fluctuations. We 
can decompose the further in trace, longitudinal and transverse-traceless 
part: 

V = -^hguv + {LC)nv + , (3.65) 

The equation to follow contains two types of covariant derivatives. Definitions 
can be found in sec. 2.1. 
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where g'^^'h^^ = ^"hj^u = and 



(LOm. = V^e. + V.e. - ^^M- ■ (3.66) 
The decomposition (3.65) is orthogonal with respect to the inner product: 
< h, h' >= J d^x^G'^^'^'h^^h'^, , (3.67) 

Here C is a constant. For positivity of (3.67) C should be greater than — |. 
Under the action of infinitesimal diffeomorphism generated by a vector the 
components of (3.65) transform as 

e^^e^ + e^, /i^/i + 2V%, hi^^h^, (3.68) 

One can fix the gauge freedom (3.68) by the condition 

^^ = 0. (3.69) 

If the background admits conformal Killing vectors (these are the vectors 
which are annihilated by the operator L (3.66)) the condition (3.69) is not 
enough and one should impose one more gauge condition on the trace part 
(see e.g. [402]). We suppose that conformal Killing vectors are absent. 

The Jacobian factor induced by the change of variables (3.65) h/^i, — > 
{h, C;., hj^^) is 

J = det(L^L)^, (3.70) 

where the determinant is calculated on the space of the vector fields (exclud- 
ing the conformal Killing vectors which we do not take into account). It is 
convenient to shift the scalar part of the metric fluctuations by V^^^ so that 
the decomposition (3.65) becomes 

V + 2V''^p)g^. + + ■ (3-71) 

Since the change of the variables h ^ a does not introduce any Jacobian 
factor we conclude that the path integral measure is 

P V = det{L^ L)^VaV^Vh^ . (3.72) 

To simplify the discussion we suppose that the background metric g^. 
satisfies the Einstein equations 

RM-^g^i'^- (3-73) 
The quadratic part of the action reads: 



29 



Due to the gauge invariance (3.74) docs not contain ^. Functional integration 
over ^ produces an infinite constant equivalent to the volume of the diffeomor- 
phism group which will be neglected. 

The kinetic term for a has a wrong (negative) sign. This represents the 
well known conformal factor problem of quantum general relativity. Different 
explanations of this phenomenon [232,375,314,348] suggest (roughly) the same 
remedy: the conformal mode must be rotated to the imaginary axis, a ia. 

The path integral can be written in terms of functional determinants 
corresponding to vector fields (indicated by the subscript V), scalars (indicated 
by S), and transverse traceless tensors (indicated by T _L): 

1 

\ — 2 

-A--j det{-Ag^pg,^ + 2Ri,p,^)--2 . (3.75) 

This expression is, in principle, suitable for calculations on some homogeneous 
spaces since the harmonic expansion on such spaces usually respects separa- 
tion of tensors to transverse and longitudinal parts. However, we must warn 
the reader again that we have neglected the presence of Killing vectors and 
of conformal Killing vectors. The ways to treat these vectors created a long 
discussion in the literature [233,128,211,8,254,396,402,335,415]. 

The representation (3.75) is not convenient on generic manifold since the 
tensor operator is restricted to the transverse modes. To remove this restric- 
tion one can either multiply (3.75) by a compensating vector determinant, 
or add a suitable gauge fixing term to the classical action (3.64) and repeat 
the quantisation procedure right from the beginning. The result can be found 
elsewhere in the literature (see, e.g., [233,128,211,396]). 

Boundary conditions for one-loop Euclidean quantum gravity must be 
diffeomorphism invariant and must lead to "good" (i.e. symmetric and elliptic) 
operators describing metric and ghost fiuctuations. The problem of finding a 
suitable set of boundary conditions in gravity appeared to be much harder 
that for the lower spin fields. There exist several proposals on the market 
[47,308,193,412,313,39,38,36,334,190]. However, neither of these proposals is 
fully satisfactory. It is clear now that proper boundary conditions must depend 
on tangential derivatives. Moreover, they probably must be non-local. 

The material of this section can be generalised to the case when an in- 
dependent torsion field is present. Explicit expressions for relevant operators 
acting on fields with different spin in the Riemann-Cartan space can be found 
in [350,132,426]. 
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4 Heat kernel expansion on manifolds without boundary 

4.1 General formulae 

For physicists, the most famihar way to calculate the heat kernel coef- 
ficients is the DcWitt iterative procedure which will be briefly described in 
sec. 4.3. Here we take a different route following the method of Gilkey. Ad- 
vantages of this method are most clearly seen on manifolds with boundaries. 
However, even if boundaries are absent we use the same method for the fol- 
lowing reasons: (i) this will ensure a smooth transition to more complicated 
material of the next section; (ii) we believe that calculations of the coefficients 
ak{f, D) arc a little bit easier even without a boundary. We are not going to 
give complete proofs of all statements. Instead, we concentrate on main ideas 
of the method. For the details an interested reader can consult the original 
paper [244] and the monographs [245,297]. 

We start with very general properties of the heat kernel coefficients. Let 
us consider a smooth compact Riemannian manifold M without boundary. To 
be able to define functions on M which carry some discrete (spin or gauge) 
indices we need a vector bundle V over M. Let D be an operator of Laplace 
type on V, and let / be a smooth function on M. There is an asymptotic 
expansion (2.21) and 

1) Coefficients with odd index k vanish, a2j+i{f, D) = 0; 

2) Coefficients a2j (/, D) are locally computable in terms of geometric invari- 
ants. 

Already the existence of the power-law asymptotic expansion (2.21) is a 
non-trivial statement. We postpone the discussion of this property to sec. 5.4. 

The second statement above is very important. It means that the heat 
kernel coefficients can be expressed as integrals of local invariants: 

ak{f,D) = trvf rx^{f{x)ak{x;D)} = Y,trv [ d^x^{fu'Ai{D)}, 

(4.1) 

where Al are all possible independent invariants of the dimension k con- 
structed from E, fl, R^upa and their derivatives. We use usual assignments 
of the dimensions when E has dimension two, any derivative has dimension 
one, etc. are some constants. For example, \i k — 2 only two independent 
invariants exist. These are E and R. Note, that we can always integrate by 
parts to remove all derivatives from /. The first statement, 02^+1 = 0, is clear 
now. One cannot construct an odd-dimension invariant on a manifold without 
boundary. 

Let us study further relations between the heat kernel coefficients which 
will turn out to generate relations between the constants . Consider now 
the case when the manifold M is a direct product of two manifolds. Mi and 
M2, with coordinates Xi and X2 respectively, and the operator D is a sum of 
two operators acting independently on Mi and M2, D — Di ® 1 -\- 1 ® D2. 
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This means that the bundle indices are also independent. As an example, one 
can consider the vector Laplacian on Mi x Af2- One can write symbolically 
exp(— tD) = exp(— tDi) ® exp{—tD2). Next we multiply both sides of this 
equation by f{xi,X2) = fi{xi)f2{x2), take the functional trace and perform 
the asymptotic expansion in t to get 

ak{x;D)= Yl ap{xi; Di)aq{x2; D2) . (4.2) 

p+q=k 

The consequences of eq. (4.2) are very far reaching. In particular, eq. 
(4.2) allows to fix the dependence of on the dimension of the manifold M. 
Consider an even more specialised case when one of the manifolds is a one- 
dimensional circle: Mi = , < Xi < 27r. Let us make the simplest choice for 
Di. Di = — All geometric invariants associated with D are defined solely 
by the D2-part. Moreover, all invariants are independent of xi. Therefore, by 
the equation (4.1), 

ak{f{x2),D)^ f <rx^Y.^Tv{f{x2)ul^Ai{D)} 

JS^xM2 J 

= 27r/ d^-'x^Y.^rv{f{x2)u\^)Ai{D2)}. (4.3) 

Here dependence of the constants uj^^-^ on the dimension n of the manifold is 
shown explicitly. On the other hand, we can use (4.2). Spectrum of the operator 
Di is known. The eigenvalues are l^, I E Z. The heat kernel asymptotics for 
Di can be easily obtained by using the Poisson summation formula 

c.W^+0(e-V*). (4.4) 



Since exponentially small terms have no effect on the heat kernel coefficients, 
the only non-zero coefficient is ao{l,Di) — y/ir. Therefore, we obtain by eq. 
(4.2) 

a,{f{x'),D) = / d^-'x^Y^Tv{f{x')ul,_i^Ai{D2)}. (4.5) 

There are no restrictions on the operator D2 or on the manifold M2. By com- 
paring the equations (4.3) and (4.5) we obtain 



= ^/4^«(n+l)• (4.6) 

This proves that the constants depend on the dimension n only through 
the overall normalisation factor (47r)^"/^ . 

To calculate the heat kernel coefficients we need also the following varia- 
tional equations (see [245], Lemma 4.1.15): 
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^|,=oafc(l, e-'^^D) = {n- k)ak{f, D) , (4.7) 
jJ,=oak{l:D-eF)^ak-2{F:D), (4.8) 
^Je=oan-2{e-''^F,e-'^fD) = 0, (4.9) 

where / and F are some smooth functions. 

To prove the first property (4.7) we note that 

4|,=oTr(exp(-e-^'-^iL>)) = Tr{2 ftD ex.p{-tD)) = -2t^Tt{f ex.p{-tD)) 
ae at 

and expand both sides of this equation in the power series in t. Eq. (4.8) can 
be checked in a similar way. To prove (4.9) consider the operator 

D{e,5)^e-^'^{D-SF). (4.10) 

We use first (4.7) with A; = n to show: 

0^jJ,=oanil,D{e,5)). (4.11) 

Then we vary the equation above with respect to 5: 

= ^|5=o^|e=oan(l,^(e,5)) = ^|.=o^|5=oa„(l,i5(6,(5)). (4.12) 

Finally eq. (4.8) yields (4.9). 

Eq. (4.8) restricts dependence of the heat kernel coefficients on the "poten- 
tial" E while (4.7) and (4.9) describe properties of the heat kernel coefficients 
under local scale transformations. 

To calculate the heat kernel coefficients we adopt the following strategy. 
First we write down a general expression for containing all invariants A^^ of 
dimension k with arbitrary coefficients . The constants are then calculated 
by using the properties derived above. The first three coefficients read [244, 
245]: 

ao(/,i^) = (47r)-"/2 / rx^trv{aof} (4.13) 

JM 

a^if, D) = (47r)-"/2l j d^x^tTv{f{aiE + a2R)}, (4.14) 
Jm 

a^if, D) = (47r)-'^/2-i- / d'^x^tvvifia^E.kk + a^RE + a^E^ 

ODD JM 

+aQR.kk + a^R? + asRijRij + a^RijkiRijki + Q;ioilijJlij)}.(4.15) 

Instead of the we use rescaled constants a/. By (4.6) the coefficients ai 
are true constants, i.e. they do not depend on n. One can check that indeed 
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no more invariants exist. For example, Rij-ij is proportional to R-jj due to the 
Bianchi identity. 

The coefficient ao follows immediately from the heat kernel expansion for 
the "free" scalar Laplacian on (see eq. (4.4)). We obtain ckq = 1- Let us 
now use (4.8). First take k — 2. Then 



If rx^trv {a,F} = [ rx^trv{F} . (4.16) 
D JM Jm 



This gives ai — 6. Take A; = 4 to see 



360 



J (Tx^tYviaiFR + 2a5FE} = ^J (txy/gtrviaiFE + a2FR} . 

(4.17) 

The equation (4.17) shows that = 180, a4 = 60q;2. 

To proceed further we need local scale transformations defined in (4.7) 
and (4.9). These scale transformations look similar to the local Weyl trans- 
formations but are not exactly the same. The scale transformations in (4.7) 
are designed in such a way that the operator D always transforms covariantly. 
This is not the case of the Weyl transformations of an arbitrary operator of 
Laplace type. For example, the scalar Laplacian (3.4) is conformally covariant 
for a special value (3.5) of the conformal coupling only. Therefore, some of 
the basic quantities are transformed in a somewhat unusual way. The metrics 
transforms as Qf^u e^^^ g^iv thus defining standard conformal properties of the 
Riemann tensor, the Ricci tensor and the scalar curvature. The functions a" 
and h in eq. (2.1) transform homogeneously. Transformation properties of uj^ 
and E are then defined through (2.3) and (2.4). One can obtain the following 
relations [244,245] 



d 

de 
d_ 

de 
d_ 

de 
d_ 

Je 
d_ 

de 
d_ 

de 
d_ 

de 



e=oRijkl — —'^fRijkl + ^jlf;ik + ^■,ikf;jl — Silf;jk — Sjkf;il , 



^^,E=-2fE+-{n-2)f,ii, 



,=oi? = -2/i?-2(n-l)/;. 



e=oE-kk = -4/-E;jfcfe - 2f.kkE + -{n — 2)f.iijj + {n — 6)f;kE.k , 
,=oRE = -AfRE + l{n- 2)f,uR - 2(n - l)f,uE , 
,=o^2^-4/£;2 + (n-2)/.,,£;. 
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d_ 

de 
d_ 

de 

d 

de 
d_ 

de 
d 

de 



e=oR;kk = -4:fR-kk - '2f;kkR " 2(n - l)f-iijj + {u - Q)f;iR;i , 

,=oR'^-AfR'-A{n-l)f.,uR, 

e=oRijRij = —4:fRijRij — 2f.iiR — 2{n — 2)f.ijRij , 

e=oRijklRijkl — RijklRijkl ^f,ijRij i 

e=o^' = -4/0' (4.18) 



Let us remind that the indices i,j,k,l are flat, so we can put them aU down 
and sum up over the repeated indices by contracting them with the Kronecker 
S. 

Let us apply (4.9) to n = 4. 



d 
de 



e=oa2(e-'^^F,e-2^^D) = 



By collecting the terms with try Ji^ d^x^{F f.^j} we obtain cci = 60:2 and 
consequently 0:2 = 1 and 04 = 60a;2 = 60. This completes calculation of 02. 

Let M = Ml X M2 with a product metric and let D = (-Ai) + (-A2) 
where Ai^2 are scalar Laplacians on Mi and M2 respectively . Eq. (4.2) yields 

04(1, -Ai - A2) = 04(1, -Ai)ao(l, -A2) + 02(1, -Ai)a2(l, -A2) + 

+ao(l,-Ai)a4(l,-A2) (4.19) 

It is clear from (3.6) that E — and Q = 0. By collecting the terms with 
R1R2 (where Ri and R2 are scalar curvatures on Mi and M2 respectively) we 
obtain 

2 /osV 



360°"^ V 6 y 

Consequently, 0:7 = 5. 

Let us apply (4.9) to n = 6. We obtain with the help of the variational 
equations (4.18) 

= try / d''x^{F{{-2a3 - lOa^ + 4a5)/.fefeE 
Jm 

+ (2a3 - 10aQ)f.iijj 

+ (2^4 - 2a6 - 20a7 - 2a8)f;nR 

+ {-8as-8ag)f.,jRij)}. (4.20) 

The coefficients in front of independent invariants in (4.20) must be zero. We 
determine = 60, ctg = 12, ag = —2 and ag = 2. 



More precisely, we assume that in (3.4) the potential [7 = 0, the conformal cou- 
pling is minimal ^ = 0, and there is no gauge coupling = on both Mi and 
M2. 
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The most elegant way to calculate the remaining constant aio is based on 
the Gauss-Bonnet theorem (see [245]). We use here a more lengthy way [337] 
which however works perfectly on flat manifolds without boundary and of 
trivial topology of M". Note that M" is non-compact. To make the heat kernel 
well defined we should suppose certain fall-off conditions on the background 
fields and on the smearing function /. The basis in the space of the square 
integrable functions is given by the plane waves exp{ikx). Therefore, for M — 
M" with flat metric the heat kernel reads: 



K{f;t)^T:tL2{fexp{-tD)) 

- jd'^^j ^e-'=nrv{/(x)exp(-iL>)e^'=n 

= / rf"^ / (^tr^{/(^) exp(t((V'' + ik^f + E))} . (4.21) 
The following integrals will be useful: 



(Pk ^_^^2 _ 1 



(27r)" (4t7r)"/2 ' 

J {2n)- (4t7r)-/2 2i^ ' 

/ -^e-'^'k^k'^k^k" = \ /o A(^^V" + g^'g"" + g^'^g"')- (4.22) 

J (27r)" (4i7r)"/2 4*2 ^ ^ ^ y y y v > 



Now we isolate exp(— t/c^) on the right hand side of (4.21) and expand the rest 
of the exponent in a power series of t. 



K{f; t)^j rx I -0^e-'^\.y \j{x) [l + ^(V^ E) - ^4(W) 
*2 

+ ^(V^V^ + V'^E + EV'^ + ^2) 

4+3 

- ^(('^V)'^; + E{kVf + {kV)E{kV)) 


4+3 

- -^((^V)'V' + V'(W)' + (W)V'(A;V)) 

+ + •••)} (4.23) 



We use the integrals (4.22) to obtain 
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^(/; ^) = J^^2 I d'-xtVy {f{x) (1 + tE 

+ ^(V V + V'^E + EV'' + 

- -(v^^; + Ev^ + v''£;v^) 

+ ^(V^V'V^V. + V'V' + VV'V^) + O(t') ) I . (4.24) 



All derivatives combine into commutators. Finally we get: 



(^^' + \E;,, + ]^^.-^"'^) + } ■ (4.25) 



Prom the equation (4.25) we conclude that aio — 30. We have obtained also 
an independent confirmation for the values of ao, cti, «3 and a^. 

The method we used above in eqs. (4.21) - (4.25) can be applied to more 
general manifolds and operators as well. The key ingredient is a convenient 
basis which should be used instead of the plane waves. For the case of a 
box with (anti-) periodic boundary condition such a basis is rather obvious 
(see [13, 14]). On curved manifold M one has to use the so-called geodesic 
waves [220,221,425,427,347] although calculations with ordinary plane waves 
are also possible [120]. 

For calculation of the coefficient we refer to [244] . The results for the 
leading heat kernel coefficients are summarised in the following equations 
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ao(/, D) = (47r)-'^/2 f d-x^tTv{f}. (4.26) 

Jm 

a2(/, i^) = (47r)--/2g-i I c?"xV^tr^{/(6E + R)}. (4.27) 

04(7, D) = (47r)-"/^360-^ / d"xV^try{/(60£;;fcfc + 60RE + 180E^ 

Jm 

+ 12-R;fcfc + 5i?^ — 2RijRij + 2RijkiRijki + SOfiyi^jj)}. (4.28) 

a6(/, £>) = (47r)-'^/2 / ci'^xV5try{^(18i?;ii,-,- + 17 R.kR.,k - ^Rij-kRij-k 
Jm il 

—4:Rjk;nRjn;k + 9Rijki,nRijkl;n + 28RR-nn — 8RjkRjk-,m 

+2ARjkRjn-kn + l^RijkiRijki-nn + 35/9i?^ — 14:/3RRijRij 
+14:/3RRijkiRijki — 208 /9RjkRjnRkn — G4:/3RijRkiRikji 

\Q j 3RjkRjnliRknli ^RijknRijlpRknlp 

— 80/9 RijknRiikpRjinp) + 360 ^ f{8ilij-k^ij-k + 2ilij-jilik;k 

-\-\2^lj-kk^ij 12^'ij^jk^ki ^Rijkn^ij^kn ^Rjk^jn^kn 

+bRVtknnkn + QE.iijj + mEE.ii + 30E.,E., + 60^3 
-\-30EVtijVtij + lOit!^/;;;;;;; + ARjkE-jk + 12R-kE-k + 30EER 
+12ER.kk + 5ER^ — 2ERijRij + 2ERijkiRijki)}- (4.29) 

Everyone who ever attempted calculations of the heat kernel coefficients 
on curved background for arbitrary spin should appreciate that the method 
presented here is a quite efficient one. One should also take into account that 
some of the universal constant were indeed calculated twice. 

The coefficients qq - are contained in [147,319]. oe was first computed 
by Gilkey [244] . The next coefficient ag has been calculated by Amsterdamski 
et al [12] for the scalar Laplacian and by Avramidi [25, 26] for the general 
operator of Laplace type. The coefficient aio has been calculated by van de 
Yen [399]. Higher heat kernel coefficients in flat space were studied in [205,206]. 

4-2 Examples 

Here we consider several simple physical systems in four dimensions and 
calculate the heat kernel coefficient 04 which defines the one-loop divergences 
in the zeta function regularization. 

4.2. 1 Yang-Mills theory in flat space 

Our first example is pure Yang-Mills theory in flat space. We are inter- 
ested in the "total" heat kernel coefficient 04°* defined by (3.53). Let us start 
with the first term describing contribution from the vector fields. We choose 
the gauge (3.44). The only non- vanishing invariants are E and Vt. The coef- 
ficient 04 is quadratic in these quantities. By using eqs. (3.46) and (3.47) we 
obtain: 
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tr^ (e^) = E^'^E^^^ = ^F^FI^Ks, , (4.30) 
trv. (a,-Q,,) = -AFlF^^Ks, , (4.31) 

where 

Ks, = ci,clp (4.32) 

is the Kilhng form of the gauge algebra. 

Let us remind that the indices k in (4.28) refer to a local orthonormal 
frame. In flat space they may be identified with the vector indices /i, u, p. In the 
case of vector fields the trace in (4.30) and (4.31) is taken over pairs consisting 
of a gauge index a, /5 or 7 and of a vector index /i, i' or p. 

For the ghost operator we have to put R = U = in (3.6) and identify the 
connection with the background gauge field according to (3.49). Consequently, 

= and 

tr^ (n!f Olf ) = -FlF^^Ks, . (4.33) 
Next we substitute (4.30), (4.31) and (4.33) in (4.28) to obtain 

ar 1 = af- 1 - 2af 1 = ^ d'x ^ F^F^^Ks, . (4.34) 

The Yang-Mills gauge group Q has usually a direct product structure, 
Q — ^1 X ^2 X ■ ■ ■ , so that on each of the irreducible components Qi the 
Killing form K is proportional to the unit matrix. Therefore, the one-loop 
divergence (4.34) reproduces the structure of the classical action with different 
charges for each Qi. We also recover the coefficient 11/3 which is familiar from 
computations of the Yang-Mills beta functions. 

4-2.2 Free fields in curved space 

Consider free quantum fields on a Riemannian manifold without bound- 
aries. "Free" means that we neglect all interactions except for the one with 
the background geometry. The heat kernel coefficients can be expressed then 
in terms of local invariants of the metric. In particular. 



28 807r2 

where C^,/po- is the Weyl tensor. 



, (4.35) 



C^.paC^'-''' = R^.paR^"'" - 2R^,R^'' + -R^ . (4.36) 

a, b, c and d are some constants depending on the spin. The first two structures 
in (4.35) which appear together with a and b are conformally invariant in four 
dimensions. This explains our choice of the basis in the space of invariants. 

The constants a, b, c and d can be evaluated by substituting particular 
expressions for E and Q obtained in sec. 3 in (4.28). Alternatively, one can use 
the analysis of Christensen and Duff [127] who calculated 04 for arbitrary spin 



39 



Table 1 

04 for various spins 



Spin 


a 


b 


c 


d 





1 


1 


30^ - 6 


m - 1/6)' 


1/2 


-7/2 


-11 


6 





1 


-13 


62 


18 





2 


212 








717/4 



fields (see also [69]). The results are collected in Table 1. Some comments are in 
order. Spin 1/2 means 4-component Dirac spinors. For spin 1 and spin 2 fields 
we took into account contributions from corresponding ghost fields. Note, that 
vector Yang-Mills fields and vector ghosts for gravity interact differently with 
the background geometry. The cosmological constant is taken to be zero. 

As a physical application of the heat kernel expansion in curved space we 
may mention, for example, the asymptotic conformal invariance phenomenon 
which was studied by using this technique in [103, 104,99, 108]. 

A similar analysis can be performed also in the presence of boundaries 
[333]. 



4.3 DeWitt iterative procedure 

The iterative method by DeWitt uses separation of the heat kernel with 
non- coinciding arguments into a leading part (which is non-analytic in t) and 
the power-law corrections. For a flat manifold this separation has been de- 
scribed in sec. 1 (cf. eqs. (1.12) and (1.13)). Let us consider a massless scalar 
field on a curved compact Riemannian manifold M. The DeWitt ansatz reads 
in this case: 

K{t; X, y; D) = (47ri)--/^A^/;^(x, y) exp (-^%^) S(i; x, y; D) , (4.37) 

where a{x, y) is one half the square of the length of the geodesic connecting 
X and y. In Cartesian coordinates on a flat manifold o"flat(a;, y) = \{x — y)^. 
^VVM is the so called Van Vleck-Morette determinant 

A,,„(.,,).£!±^E±Z)). (4.38) 

As a consequence of the heat equation (1.10) the kernel H should satisfy 
(dt + t-\W^a)V^ + Ay'JlDA]!}^) E = (4.39) 
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with the initial condition 



E{Q-x,y-D) = l. (4.40) 

The essence of the DeWitt method is to look for a solution of the equation 
(4.39) in the following form 

oo 

E{t-x,y-D) = Y.t^h2j{x,y-D). (4.41) 
The initial condition (4.40) yields 

6o = l- (4.42) 

The recursion relation 

(j + {V^a)V,) + A-\/lDA\!}j^b2U-i) = , (4.43) 

which follows from (4.39) and (4.41), allows, in principle, to find the coinci- 
dence limits X = y oi higher heat kernel coefficients b2j- An important ingredi- 
ent of such calculations is the coincidence limits of symmetrized derivatives of 
the geodesic interval a{x, y) [395,55]. This method becomes very cumbersome 
beyond 04. A refined nonrecursive procedure to solve the DeWitt equation 
(4.43) was used by Avramidi [26] to calculate ag (see also [34] for a short 
overview). The method of DeWitt can be naturally extended to treat coinci- 
dence limits of the derivatives of b2j{x, y; D). 

The recursion relations can be generalised to the case of manifolds with 
boundaries [315,318]. However, for a practical use the functorial methods of 
sec. 5 seem to be more convenient. 

4^.4. Non-minimal operators 

Quantisation of gauge theories oftcnly leads to second order differential 
operators which are not of the Laplace type. For example, by taking k 7^ 1 in 
(3.52) one obtains the following operator acting on the gauge field fiuctuations: 

Djrl = -Ag,, + (1 - ^) ^^^^ ' (4-44) 

where, for simplicity, we suppose that the manifold M is fiat {R^upa = 0) and 
the gauge group is abelian (c^^^ = 0). The leading symbol of the operator (4.44) 
(the part with the highest derivatives) has a non-trivial matrix structure. Such 
operators are called non-minimal. 

In some simple cases (see, e.g., [240,7,411]) the spectral problem for non- 
minimal operators can be reduced to the Laplacians. If the leading part of a 
non- minimal operator has a form similar to (4.44) but the lower order part 
is more or less arbitrary , necessary generahsations of the DeWitt technique 

This case covers most of the physical applications. 
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were suggested by Barvinsky and Vilkovisky [55]. The technique was further 
developed in [269,123]. Complete calculation of 04 required extensive use of the 
computer algebra [270] . Most general non- minimal operators were considered 
in [32] where, because of great technical complexity, only the first two heat 
kernel coefficients were analysed. Some calculations in various physical systems 
with non-minimal operators can be found in [209,357,266,361]. 
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5 Heat kernel expansion on manifolds with boundaries 

5. 1 Two particular cases 

We start our analysis of manifolds with boundaries with two simple ex- 
amples. First, let us consider a one-dimensional manifold M = [0,7r]. Let 
D = — 19^. We consider both Dirichlet (2.15) and Neumann (2.16) boundary 
conditions (taking 5" = in (2.16) to simplify the calculations). The eigen- 
functions of D are: 

Dirichlet: sin(/a;), / = 1, 2, ... ; 

Neumann: cos(/x), I — 0,1,2, . . . (5.1) 

The eigenvalues are in both cases, where I is a positive integer for Dirichlet 

boundary conditions, and I is a nonncgativc integer for Neumann boundary 
conditions. The heat kernel asymptotics can be calculated with the help of the 
Poisson summation formula (4.4): 

K{t,-dl,B-)=Y.eM-tl') = l(J^-l)+0{e-'/') , (5.2) 
i>o z \ \ I / 

K{t,-dl,B^)=j:eM-tl') = UJ^+l)+0{e-'/') , (5.3) 

where for the later use we explicitly mention the boundary operators B~ and 
which define Dirichlet (2.15) and Neumann (2.16) boundary conditions 
respectively. 

Let us modify a little bit the example above by allowing for a non-zero S 
at one of the components of the boundary: 

9.0|.=o = 0, (9, - 5)0U=, = (5.4) 

Note, that at x = tt the derivative with respect to an inward pointing unit 
normal is —d^- The eigenf unctions 

(pk = cos{kx) (5.5) 

satisfy the boundary conditions (5.4) a,t x — 0. The spectrum is defined by 
the condition at x — tt which reads : 

— A;sin(A;7r) = 5'cos(A;7r) . (5.6) 

In this example we restrict ourselves to the linear order in S. We suppose that 
S is small and that the spectrum (5.1) is only slightly perturbed: 

ki^l + zi. (5.7) 
For S > also a negative mode cosh(fcx) appears. Here we take S <0. 
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The equation (5.6) gives in this approximation 



zi — — for I > , 

zl^-- for / = 0. (5.8) 

TT 

We expand also the heat kernel in a power series in S and use (5.8) and (5.2) 
to obtain: 



1>Q 



l>0 \ TT J TT 

where we also dropped exponentially small terms in t. 

Our next example is the heat kernel expansion for scalar fields in a ball 
with Dirichlet and Neumann boundary conditions [392,416,291,328,72]. The 
metric of the unit ball in reads: 

C?s2 ^ ^^2 ^ j.2^q2^ < r < 1, (5.10) 

where dQ'^ is the metric on unit sphere S*""^. The scalar Laplace operator has 
the form 

A0 = (d^^ + ^^a, +("-^) a) , (5.11) 

where ("^i)A is the Laplace operator on 5"""^. The eigenfunctions of the op- 
erator (5.11) are well known 

<Pi,, oc rP-")/V(„_2)/2+KrA^/2)F(o(a;'^). (5.12) 

Jp are the Bessel functions. The eigenvalues —A are defined by boundary condi- 
tions. Y(;) (x*) are n-dimensional scalar spherical harmonics. Their degeneracies 
arc 

_ {2l + n-2){l + n-3)\ 

n^^\ ■ ^^-^^^ 

To proceed further it is convenient to consider the zeta function which 
may be presented through a contour integral [72] : 



°° r rik B 

C(«.^) = g^./^«=-^-a^*(.-.,/.+.(*). (8.14) 



A general discussion of the representation of the zeta function by contour integrals 
can be found in [299]. 
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where is a function which has zeros at the spectrum k — y/X. For Dirichlet 
boundary conditions this function reads: 

= k-''Mk) . (5.15) 

For Robin boundary condition is given by a somewhat more comphcated 
combination of the Bessel functions. The contour 7 runs counterclockwise and 
encloses all the solutions of $y = on the positive real axis. Note the presence 
of k~'^ in (5.15) which is included to avoid unwanted contributions coming from 
the origin k — 0. 

The next step is to rotate the contour to the imaginary axis and to calcu- 
late residues of r(s)(^(s, D) (as prescribed by (2.26)). It is interesting to note 
that the heat kernel coefficients are defined by several leading terms in the 
uniform asymptotic expansion of ^^{iuz) for large v and fixed z. For further 
details we refer to [72] . 

Euclidean ball was frequently used in calculations of the heat kernel co- 
efficients and functional determinants. Apart from the papers already quoted 
above also the computations for a scalar field [163,75,158,73,339], spinors 
[141, 288, 289, 159, 164, 298, 179, 199], abelian gauge fields [307, 196, 192, 404, 
410, 159, 298, 194, 197, 181, 195, 64] should be mentioned. A similar tech- 
nique works also for more complicated geometries as, e.g., the spherical cap 
[52,287,51,288,163]. Note that here we consider local boundary conditions 
without tangential derivatives only. For other types of boundary operators see 
sec. 5.4. 

5.2 Dirichlet and Neumann boundary conditions 

Let us now find analytic expressions for the heat kernel coefficients in 
terms of geometric invariants. Here we follow the method of Branson and 
Gilkey [90, 245] . For both Dirichlet and modified Neumann (Robin) boundary 
conditions the heat kernel coefficients are locally computable. This means that 
ttk may be represented as a sum of volume and boundary integrals of some 
local invariants: 




(5.16) 

where f^^^ denotes j-th normal derivative of the smearing function /. Note, 
that we cannot now integrate by parts to remove normal derivatives from /. 
This refiects the distributional nature of the heat kernel asymptotics. 



Some of these works consider contributions of the so called physical modes only. 
As explained in [413] complete answer for the effective action and for the scaling 
behaviour must include contributions from ghosts and non-physical modes. This 
applies also to the spherical cap case considered below. 
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The volume terms ak{x, D) are the same as in the previous section (see 
(4.26) - (4.29)). Here we evaluate the boundary terms a^jix, D, B). Obviously, 
the canonical mass dimension of akj is A; —j — 1. In addition to the usual bulk 
invariants E, R, etc, a^j can also contain specific boundary quantities as the 
extrinsic curvature Lab or S (for Robin boundary conditions). Note, that Lab 
and S arc defined on the boundary only and, therefore, can be differentiated 
only tangentially ^ . Canonical mass dimension of Lab and S is +1. 

In this section we explicitly calculate the first three heat kernel coefficients 
Qq, ai and a2- Basing on the considerations of the preceding paragraph we may 
write: 

ao(/,i^,^^) = (47r)-"/2 / rx^trvif), (5.17) 

JM 

ai{f,D,B^) = (47r)-("-i)/2 f ^-^^^bftrvif), (5.18) 

JdM 

a2{f,D,E^) = ^(47r)-"/2{^(i"xV^tr^(6/£; + /i?) 

+ / r-^xVhtrv{bifLaa + bff.,n + b^ f S)] , (5.19) 

JdM ) 

where B"^ denotes either Robin (S+) or Dirichlet {B~) boundary operator. 6=*= 
are some constants. To keep uniform notations we formally included S also in 
the expression for Dirichlet boundary conditions. It should be assumed that 
5 = for that case. 

Since the product formula (4.2) is still valid with obvious modifications 
for the boundary contributions, we can again consider the case M — Mi x 
and repeat step by step the calculations (4.3) - (4.6) thus arriving to the same 
conclusion that the constants b"^ do not depend on dimension of the manifold 
(all explicit n-dependence of the heat kernel is contained in the power of 47r). 
This property is crucial for our calculations. 

The one dimensional examples (5.2) and (5.3) immediately give 

bt^-bT^\. (5.20) 
Our next example (5.9) controls linear terms in S and gives 

bt^ 12. (5.21) 

For Dirichlet boundary conditions we have set S = 0. Therefore, 64 plays no 
role. 

To define the constants bf and bf we use the conformal variation equation 
(4.7). Variations of E and R are given by (4.18). Conformal properties of the 

^^See sec. 2.1 for more information on differential geometry of manifolds with 
boundary. 

For M = [0, tt] the boundary integral reduces to a sum of two contributions from 
X = and x = n. 
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second fundamental form Lab are dictated by that of the metric: 



de 



-0Laa^-fLaa-{n-l)f;n. (5.22) 



We wish to keep the boundary conditions invariant under the conformal trans- 
formations. This imphes that the boundary operators must transform ho- 
mogeneously. This property holds automatically for Dirichlet boundary condi- 
tions (2.15). In the case of modified Neumann boundary conditions conformal 
transformation of S should cancel inhomogeneous term in the connection a;„. 
This yields 

^J,=oS = -fS + ^in-2)f,n. (5.23) 

Next we substitute the variational formulae (4.18), (5.22) and (5.23) in (4.7) 
with k = 2. We collect the terms with /.„ on the boundary to obtain 

-{n - 4) - (n - 1)62" = (n - 2)63 (5.24) 

for Dirichlet conditions and 

-{n - 4) - (n - l)bt + 6(n - 2) = (n - 2)6+ (5.25) 

for generahsed Neumann (Robin) ones. Since the constants 6+ do not depend 
on n the two equations (5.24) and (5.25) are enough to define these constants: 

6+ =62 =2, 6+ = -63- = 3. (5.26) 

This completes the calculation of 02. 

The coefficients 03 and 04 can be obtained as particular cases of more 
general formulae of the next subsection. 



5.3 Mixed boundary conditions 



Let us now turn to mixed boundary conditions (2.17) which, as we have 
seen in sections 3.3 and 3.4, are natural boundary conditions for spinor and 
vector fields. These l)oundary conditions depend on two complementary local 
projectors n_ and 11+ = 1 — n_ which define subsets of components of the 
field satisfying Dirichlet and Robin boundary conditions respectively . More 
precisely, 

n_0|9M = O, (Vn + 5)n+0|aM = O. (5.27) 

Consequently, there is one more independent entity (as compared to Robin 
case, for example) on which the heat kernel coefficients for mixed boundary 
conditions can depend. This makes the calculations somewhat more compli- 
cated. It is convenient to define 

X = n+ - n_ . (5.28) 

Examples of such projectors are given by eq. (3.35) and below eq. (3.58). 
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Calculation of the coefficients up to A; = 4 can be found in [90] (see 
also [404] for some corrections). The result reads: 



ao{f,D,B) = (47r)-"/^ / d^x^trvif). (5.29) 
ai(/,L',B) = i(47r)-("-^)/2 / r-'x^/htrvixf). (5.30) 

4 JdM 

+ [ (T-^xy/htiviVLaa + 3x/n + 12/5)1 . (5.31) 

JdM ) 

as{f,D,B) = ^(A7r)-^''-'y'l^^r-'xVhtrv{f{96xE+l&xR 

Lab Lab 

+96SLaa + 192^2 - 12X:aX:a) + /;n((6n+ + 30n_)L„„ 
+96S)+2Axf;nn}. (5.32) 

a4(/, D, B) = -i-(47,)-V2| j rx^trv{f{60E,ii + GORE + 180E^ 
obU Jm 

+30QijQij + 12R.ii + 5R — 2RijRij + 2RijkiRijki)} 
+ J ci"-^a;\//itry{/{(240n+ - 120n_)£;.„ 

+ (42n+ - 18n_)i?;„ + 24L,,;fef, + OLab;ab + 120ELaa 
+20RLaa + ARananLbb — ^'^RanbnLab + ^RabcbLac 

+^{(280n+ + 40n_)L„„L66Lee + (168n+ 

-2UIl^)LabLabLcc + (224n+ + 320n_)L„6L6eL„c} 
+7205^; + 1205i? + OSRanan + UASL^U^ + A%SL,j,L^ 

+mS''Laa + 480^3 ^ ^2QS.,aa + Q^XX An ' l'^X:aX:aLbb 
-24X:aX:bLab - 120X:aX:aS} + f-nC^SOxE + 30xR + ORanan 

+^{(84n+ - mu_)LaaLbb + (84n+ + 6on_)L„bLa&} 

+72SLaa + 240-5" - 18X:aX:a) + /;nn(24L,„ + 120-S) 
+30x/;nn}}. (5.33) 

Dirichlet and modified Neumann (Robin) boundary conditions are recovered 
when n+ = or n„ = respectively. 

The coefficients with k = 0,1,2,3,4 were calculated in [3f 9, 292, 388, 
320,329,90,315,318,404]. A different algorithm was suggested in [131]. For 
"pure" (not mixed) boundary conditions the coefficient 05 was calculated by 
Branson, Gilkey and Vassilevich [93] for the special case of a domain in flat 
space or of a curved domain with totally geodesic boundaries. Kirsten [296] 
generalised these results for arbitrary manifolds and boundaries. Branson, 
Gilkey, Kirsten and Vassilevich [91] calculated the coefficient for mixed 
boundary conditions. 
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5-4 Other boundary conditions 

From the technical point of view boundary conditions for a Laplace type 
operator are needed to exclude infinite number of negative and zero modes 
and to ensure self-adjontness. In principle, any linear relation between the 
boundary data (f)\dM and (j);n\dM is admissible as long as it serves this purpose. 
Here we consider two physically motivated examples of boundary conditions 
which contain tangential derivatives on the boundary. These examples should 
give the reader an idea of what can be expected for a more general boundary 
value problem. 

5.4-1 Boundary conditions with tangential derivatives and Born-Infeld ac- 
tion from open strings 
The boundary condition 



is the simplest condition containing both normal and tangential derivatives. 
Va and S are some matrix valued functions defined on the boundary. Such 
or similar structures appear in open strings (cf. sec. 3.2) and in chiral bags 
(cf. sec. 3.3). They also describe photons with the Chern-Simons interaction 
term concentrated on the boundary [184, 83] and may be relevant for solid 
state physics applications. The conditions (5.34) appeared in the mathematical 
literature [257,243]. 

Several heat kernel coefficients for the boundary conditions (5.34) (which 
are called oblique) have been calculated by McAvity and Osborn [316] and by 
Dowker and Kirsten [167,168]. Avramidi and Esposito [37,36,38,35] lifted some 
commutativity assumptions and proved a simple criterion of strong ellipticity 
(see below). 

Consider a differential operator Q, not necessarily of the second order. Let 
is separate the part containing highest derivatives and replace 9^ ik^, like 
in doing the Fourier transformation. In this way we obtain an object Aq = 
a^'^^^■p{x)k^ky . . . kp which is called the leading symbol of Q. If the operator 
Q acts on fields with indices (i.e. in a vector bundle), the index structure of 
Q is inherited by Aq. In other words, for fixed k the leading symbol A is 
a matrix valued function (an endomorphism of the vector bundle). If Aq(x) 
is non-degenerate for all A; ^ 0, the operator Q is called elliptic. For the 
Dirac operator Q = p the leading symbol is Aq = '^^k^. Therefore, p is 
obviously elliptic. The Laplace type operators are also elliptic since A^ — k"^. 
Ellipticity means that "at large momenta" the operator Q is dominated by 
its' highest derivative part. This is the property which guarantees that on a 
compact manifold without boundaries Laplace operators have at most finite 
number of negative and zero eigenvalues and which ensures existence of the 
heat kernel. On manifolds with boundaries just ellipticity is not enough to 




(5.34) 



49 



ensure nice properties of the spectrum. There is an additional requirement, 
called strong ellipticity, which should be satisfied by the boundary operator 
(see [245] for details). Dirichlet and Robin boundary conditions are always 
strongly elliptic. Oblique boundary conditions are strongly elliptic if and only 
if |r^| < 1. If this inequality is violated, infinite number of negative eigenmodes 
appears (a simple example can be found in Appendix B of [306]). Note that 
the boundary conditions (5.34) correspond to a symmetric D if the matrices 
Ta are anti-hermitian. Therefore, = FaF" is typically negative. 

Another complication stems from the fact that F^ is dimensionless. Con- 
sequently, arbitrary powers of F^ may enter Ofc, so that instead of the unde- 
termined constants (cf. (5.17) - (5.19)) one has to deal with undetermined 
functions of F. The problem becomes more tractable if we suppose that F^'s 
commute among themselves, [F^, F;,] = 0. For this case, the coefficients Oi and 
02 have been calculated by McAvity and Osborn [316] (oq is still given by 
(5.17)): 



(5.35) 



a^{f, D) = (47r)-("-i)/2 / d^-'x^iYy{f^{V)) , 

JdM 

a^if, D) = ^(47r)-'^/2 { / (Tx^iivi&fE + fR) 
+ / (r-^xy/hiYv{ho{V)fLaa + 6l(F)/;„ 

JdM 



(5.36) 



where 



^=4 



ho 



1 



- 1 



1 + F2 
6 



i=artanh (v'^) 



hi — ^ "^^ artanh ^V— F^j — 3 , 
62 



+ 2, 



(5.37) 



12 



1 + F2 ' 
<7=^(2-M- 

We see, that if F^ approaches —1 the heat kernel blows up indicating violation 
of the strong ellipticity condition. 

As an example, let us consider the open string sigma model of sec. 3.2. 
To simplify the subsequent calculations we put Bab = 0. In the zeta function 
regularization divergent part of the effective action (2.31) reads (l/s)a2(-D) = 
(1/5)02(1, -D). On general grounds we expect that 02 repeats the structure of 
the classical action (3.10): 



02 



dM 



dTp^fdrX^ + bulk terms. 



(5.38) 
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where is a beta function ^° . We put Lab — (otherwise we have had to 
introduce a dilaton coupUng on the boundary to achieve renormahzabihty) . We 
also suppose that the target space metric Gab is trivial. With these simplifying 
assumptions the beta function can be easily calculated from (5.38), (5.36), 
(5.37), (3.21): 

4^' = -^(5aFbc)(1 + F^)ba- (5.39) 

By a lengthy but straightforward calculation one can demonstrate that the 
condition /?^^ = is equivalent to the equations of motion following from the 
Born-Infeld action on the target space: 

£bi = J dX^det{l + iF) = J dXexp Qtrln(l + F^)^ , (5.40) 

where i appeared due to our rule of the Euclidean rotation for the gauge fields. 

The Born-Infeld action has been derived from the beta functions of the 
open string sigma model in [1,112] confirming an earlier work [212] which used 
different methods. The heat kernel analysis was performed by Osborn [352] 
and then repeated in [306] for more general couplings. 



5.4-2 Spectral or Atiyah-Patodi-Singer (APS) boundary conditions 

Spectral boundary condition were introduced by Atiyah, Patodi and Singer 
in their study of the Index Theorem [19,20,21]. These boundary conditions 
are global, i.e. they cannot be defined by using local data only. 
Consider a Dirac type operator 

P = ^7''V/, + S , (5.41) 

where V is a covariant derivative with a compatible connection: 

V^7. = 0, (5.42) 

i.e., the gamma matrices are covariantly constant, f is a zeroth order opera- 
tor (a matrix valued function). We also suppose that the connection in V is 
unitary. This means that the connection one-form is represented by an anti- 
hermitian matrix in a suitable basis. We restrict ourselves to the case when 
= £, so that the operator^ is formally self-adjoint in the bulk. Note, that 
compatible unitary connection is not unique. Consider a first order differential 
operator on the boundary: 

P = 7n7"V„ + {£-in - inS) + e(x) , (5.43) 

In this case the subscript A \s a target space vector index, while the superscript 
[A\ indicates the coupling. Here we do not consider other beta functions Pi^^\ /J'^^ 
etc. 
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where 0(,x) is a hermitian matrix valued function on dM. The operator P is a 
self-adjoint operator of Dirac type on the boundary. All functions on dM can 
be decomposed in positive, negative, and zero modes of P. Let us define n_ 
as a projector on the space spanned by non-negative eigenspaces of P. Then 
the equation 

n_0|aM = (5.44) 

defines the APS boundary conditions. 

Of course, relative complexity and non-locality limits physical applica- 
tions of spectral boundary conditions. However, they appeared in a number 
of axial anomaly calculations [369,280,278,342,341,312,208,201], quantum 
cosmology [142], and in works on the Aharonov-Bohm effect [327,61]. In brane 
models spectral boundary conditions describe T-selfdual configurations which 
may be interpreted as mixtures of D-branes and open strings [406]. 

General form of the heat kernel expansion for spectral boundary condi- 
tions was established by Grubb and Seeley [260,261,262]: 

n— 1 oo 
k=0 j=n 

In the contrast to all previous cases the expansion (5.45) contains log- 
arithms of the proper time t. Although, such terms appear "typically" [241] 
for pseudo-differential operators, they may lead to rather unpleasant physical 
consequences. As follows from (2.25), non-zero means that the zeta function 
has a pole at s = 0, and, therefore, the expression (2.32) for the renormalised 
effective action does not make sense. Fortunately, for the APS boundary con- 
ditions = if / = 1 near the boundary. Consequently, the integrated zeta 
function is regular at s = and eq. (2.32) still can be used. However, all calcu- 
lations involving localised heat kernel coefficients remain problematic. Many 
logarithmic terms vanish if the manifold M has a product structure near dM. 
In non-product cases strong criteria of "partial vanishing of logarithms" have 
been found recently [263]. 

Another problem with the heat kernel expansion (5.45) is that the coef- 
ficients may have a more complicated dependence on n than just a power 
of (47r). After this long list of troubles it is not a surprise that only the coef- 
ficients Ofe with k < n are locally computable. The coefficient oq is given by 
(4.26). Oi and 02 have been calculated in [172]: 



Let us remind that to prove the simple dependence on n we used the product 
formulae (4.3), (4.4). We assumed that the spectral problem can be "trivialised" 
in one direction. In the present case all tangential coordinates enter the operator 
(5.43) on equal footings. Therefore, the proof does not go through. 
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ai = (47r)(i-")/2l(/3(n 



1) / (T-^xVhXYvif) , 

JdM 



a2 = (47r)-"/2 



n 



2{n - 2) 



1 - ^^P{n) ) /;n 



where 



n 



-1 



(5.46) 



(5.47) 



The coefficient can be found in [242]. 

Some string theory apphcations suggest [406,414] that spectral boundary 
conditions can be defined directly for a second order differential operator. Ex- 
istence of the asymptotic expansion (5.45) and vanishing of leading logarithms 
for such problems have been stated in [264]. 
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6 Manifolds with singularities 

All results on the heat kernel expansion formulated in the previous sec- 
tions are valid on smooth manifolds only. If there are boundaries, they also 
have to be smooth. As well, any singularities in the potential term or in the 
field strength are strictly speaking forbidden. However, many physical models 
deal with singular backgrounds. Even if such backgrounds may be represented 
through certain limiting procedures from smooth configurations, the heat ker- 
nel coefficients are not given by limits of their "smooth" values. The most 
visible manifestation of failure of the smooth field approximation is that 
with sufficiently large k arc divergent. Usually, the presence of singularities 
changes even the structure of the heat kernel kernel expansion as compared 
to the smooth case. 

6.1 Non-integrable potentials 

According to (1.21) divergences in the effective action are defined by 
integrated heat kernel coefficients. Although the formulae (4.26) - (4.29) for 
the localised heat kernel coefficients are valid on non-compact manifolds (pro- 
vided the smearing function / falls off sufficiently fast), transition to the 
integrated heat kernel is not that straightforward. Already the coefficient 
ao(-D) = 00(1, D). which is proportional to the volume, is divergent. This 
divergence is usually removed by replacing det{D) in (1.6) by 

det(D)/det(Do), (6.1) 

where the operator Dq = —d^ + describes a free particle propagation in 
an "empty" space. It is argued that since Dq does not depend on "essential" 
variables division by det(Do) does not change physical predictions of the the- 
ory. In all subsequent formulae the heat kernel K{t; x, y; D) is then replaced 
by the subtracted heat kernel 

K,nh{t; X, y) = K{t; x, y; D) - K{t; x, y; Do) . (6.2) 

In fiat space the coefficient Oq corresponding to Kguh is identically zero . If 
the field strength and the (subtracted) potential E + have a compact 
support or decay sufficiently fast at the infinity, the small t asymptotic expan- 
sion of Ksuh{t', X, x) is integrable on the whole M. If not, the very structure of 
the global heat kernel may be changed. 

As an example of non-integrable potentials consider the harmonic oscil- 
lator in one dimension. The Schrodinger operator reads 

D^-dl + u'^x^ . (6.3) 

On a curved manifold the subtraction procedure is more subtle. On has to define 
a reference metric which differs from the physical one on a compact submanifold. 
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If we consider the problem on the whole real axis, M = R, the potential term 
is not intcgrablc. Already the expression (4.27) for 02(1, -D) diverges. There- 
fore, analytic expressions of sec. 4 cannot be used in this case. However, the 
(integrated) heat kernel can be easily calculated. Eigenvalues of the operator 
(6.3) are contained in almost any textbook on quantum mechanics: 

A,- = i/(2j + 1), j = 0,l,2,... (6.4) 

The integrated heat kernel reads 

00 -1 
K{t; D) = Y^ e-*"(2^+^) = ^[smh{ut)]-^ . (6.5) 

As i — > it behaves like 1/t while for smooth rapidly decaying potentials in 
one dimension the leading singularity in the heat kernel is 1/ Vi. This state- 
ment may be generalised to higher dimensions, li D = —d^ + P^^x^x^ with 

a non-degenerate matrix P^j^i, on M — M", the leading term in K{t; D) is 
(2t)-'^(detP)-^/2 [29]. 

6.2 Conical singularities 

Conical space is defined as M = [0, 1] x A^ where A^ is an n— 1-dimensional 
manifold called the base. The metric of the cone has the form 

{dsf = dr'' + r'^dn'' , (6.6) 

where r G [0, 1] and dQ"^ is the line element on the base A^. This metric is, 
in general, singular at r = 0. However, if we take the unit sphere S"'^^ = N 
with standard round metric, the singularity disappears and we obtain the 
n-dimensional unit ball (5.10). If a manifold has singular points where the 
metric can be approximated by (6.6) we say that this manifold has conical 
singularities. 

Conical singularities appear in many physical applications. First of all, 
with N = the metric (6.6) is the Euclidean version of the Rindler metric. 
Conical singularities appear in classical solutions of the Einstein equations 
[389,217,45] and in the supermembrane theory [400]. Gravitational field of a 
point mass in three dimensional gravity is a conical space [146,145]. There 
arc evidences [376] that "conifolds" dominate the path integral for quantum 
gravity in topological sectors. 

Sommerfeld [391] was probably the first to consider the heat kernel in 
the presence of conical singularities. The mathematical theory of the heat 
kernel asymptotics with conical singularities was developed almost 100 years 
later [122,97,98]. There two peculiar features of these asymptotics. First, 
the heat kernel expansion contains in general both integer and half-integer 
powers of t even without boundaries. Second, a non-standard Int term may 
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be contained in the asymptotic series . 

On a manifold with conical singularities no closed analytical expression 
for the heat kernel coefficients is available. However, usually it is possible 
to disentangle contributions of the singularities from the smooth part. For 
example, ii N — S^, dD,^ — d(fi^ with (f e [0, a] only 02 receives a contribution 
from the tip of the cone: 

47r2 - q;2 

In many particular cases of conical singularities a very detailed analysis 
of the heat kernel expansion has been performed [154,155,121,156,226,227, 
129, 130, 230, 85, 140, 385, 169, 66]. One-loop computations on general orbifolds 
were considered recently in [255]. 

6.3 Domain walls and brane world 

Delta function is an example of an extremely sharp background potential. 
Let us consider a manifold M and a submanifold E of the dimension n — 1. 
Let 

D[v] ^D + v5^. (6.8) 

D is an operator of Laplace type (2.2). Let h be the determinant of the induced 
metric on E. Then 5-s is a delta function defined such that 

/ dx^6j:f{x) = [ dxVhfix) . (6.9) 

The spectral problem for D[v] on M as it stands is ill-defined owing to 
the discontinuities (or singularities) on E. It should be replaced by a pair of 
spectral problems on the two sides of E together with suitable matching 
conditions on E. In order to find such matching conditions, we consider an 
eigenfunction (px of the operator (6.8): 

D[v](l>x^X(l>x- (6.10) 

Let us choose the coordinate system such that e„ is a unit normal to E and 
= on E. It is clear that 0;^ must be continuous on E: 



|a:"=+0 



(6.11) 



Otherwise, the second normal derivative of (px would create a S' singularity on 
E which is absent on the right hand side of (6.10). Let us integrate (6.10) over 
a small cylinder C = C"~^ x [— e, -|-e] 

/ d^x^ (-V^0A - [vl<Px + {E + X)<Px]) + f (r-^xVhv<t>x = . (6.12) 



The present author is not aware of any simple example where the In t terms 
actually appear. 
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We now take the limit as e ^ 0. Since the expression in the square brackets 
in (6.12) is bounded, the contribution that this term makes vanishes in the 
hmit. We obtain 

0=/ (^-^x^/h{-Wn(^)x\xn=+o + Vn(^)x\xn=-o + v(^)x) ■ (6.13) 

J c 

Since C and A are arbitrary, we conclude that a proper matching condition 
for the normal derivatives is 

-V„0|a;n=+o + V„0|a;n=_o + V(p^O. (6.14) 

Physically this problem corresponds to two domains separated by a pen- 
etrable membrane S (a domain wall). In many cases penetrable membranes 
are better models of physical boundaries then just boundary conditions which 
are imposed on each side of E independently and thus exclude any interac- 
tion between the domains [285,2]. 5-potentials are being used in quantum 
mechanics [6] where one studies the Schrodinger equation (which is nothing 
else than the imaginary time heat equation). The Casimir energy calculations 
have been performed e.g. in [374]. In the formal limit v ^ oo one obtains 
Dirichlet boundary conditions on E, although the heat kernel coefficients are 
divergent in this limit (see below). 

Further generalisations are suggested by the brane world scenario [363, 
362] which assumes that our world is a four dimensional membrane in a five 
dimensional space . According to the Israel junction condition [283] the met- 
ric in such models cannot be smooth on E. Typical form of the metric near E 
is 

{dsf = {dx'^f + e-"l^"l(^is„_l)^ (6.15) 

where a is a constant and where (dsn-i) is a line element on the (n — 1)- 
dimensional hypersurface S. Due to the presence of the absolute value of the 
n-th coordinate in (6.15), the normal derivative of the metric jumps on S. One 
can think of two smooth manifolds M+ and M~ glued together along their 
common boundary E. Neither Riemann tensor, nor matrix potential E must 
be continuous on S. Also, the extrinsic curvatures L^^ and of E considered 
as a submanifold in and in M~ respectively are, in general, different. All 
geometric quantities referring to M~ (respectively, M") and their limiting 
values on E will be supplied by a superscript (respectively, "— "). 

For the case at hand there is still an asymptotic expansion (2.21) for the 
heat kernel. The heat kernel coefficients can be decomposed as 

akif, D[v]) = at if, D) + a^(/, D)- + a^(/, D, v) , (6.16) 

where at{f,D) are known volume contributions corresponding to (cf. 
(4.26) - (4.29)). The coefficients are given by integrals over E of some local 



A similar scenario was proposed earlier in [371], sec [370] for a review. 
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invariants. Note, that 



^a^Vt- V- (6.17) 

being a difference of two connection is a (pseudo-) vector with respect to all 
space-time and gauge symmetries. Therefore, uja can be used for constructing 
the surface invariants. 

To make the formulae more symmetric we introduce two inward pointing 
unit normals v'^ and v~ to E in M+ and M~ respectively. We do not suppose 
that the smearing function / is smooth on E (therefore, there is no relation 
between /.,^+ and but we assume continuity of f-f^ = f~ = f on E. 

The surface invariants can be constructed from L^^, Rfji^i, -E^, v, Ua and 
their derivatives. This gives much more invariants than wc usually have for 
a boundary value problem. There are, however, some properties of which 
simplify the calculations considerably. First of all, must be invariant with 
respect to interchanging the roles of M"*" and M~. Also, must vanish when 
the singularity disappears. The first property excludes, for example, the term 
f{E^ — E~) which changes sign under M"*" ^ M~ . The second requirement 
excludes the invariant f{E^ + E~) because it survives even if there is no 
singularity on E. These simple arguments show that af does not contain E^ 
even though such terms are allowed on dimensional grounds. It is also very 
helpful that in some particular cases the problem in question can be reduces 
to a sum of Dirichlet and Robin boundary value problems [82,248]. 

The coefficients a|^, /c = 0, 1, 2, 3, read 

a^(/,A^) = 0. 
af(/,L',^) = 0. 

(/,D,^) = (47r)-"/2l Jj--'xVhtTy{2f{Ll + L'J - 6fv}. (6.18) 
af (/, D, v) = (47r)(^-'^)/^^ r-'xVhtry {^/(L+ L+ + L^L;, 

ab) 

+9(L+ + L- )(/J+ + /-_) + 48/^;^ + 24/a;„a;„ 
-24/(L+ + L-Jv - 24(/++ + r^.)v] . 

The coefficients and are too long to be presented here in full generality. 

Therefore, we restrict ourselves to the case of smooth geometry {R^^i = Rijkii 
L'l^ = —L~^), smooth connection {oJa = 0), and smooth smearing function 
(/;!,+ = = /;n)- In other words, the only singularity comes from the 

surface potential v. 
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af (/, D, v) = (47r)-'^/2 ^n-i^^^^^ |_ 1 ^^3 _ 1^^^ _ 

af (/, D) = (47r)-('^-^)/2 d'^'^xx/^trv^ V + ^fRv' (6.19) 

1 5 5 2 1 2I 

The heat kernel coefficients for E = 5"""^ C and = const, were calcu- 
lated in [77]. Generic E with arbitrary v was considered in [82]. Moss [330] 
added a non-smooth connection (sec also [173]). Calculations on a particular 
brane-world background can be found in [343]. The heat kernel coefficients in 
the general setting described here were calculated in [248]. This latter paper 
also considered renormalization of the brane-world scenario and predicted a 
non-standard Higgs potential on the brane. Related calculations in wormhole 
models were done in [294,293]. The C function for brane-world geometries with 
matching conditions (6.11), (6.14) was considered recently in [186,332] (where 
one can find some further references). 

It is very well known [6] that the conditions (6.11), (6.14) (which we 
call transmittal) are not the most general matching conditions which can be 
defined on a surface. In general, boundary values of a function and of its' 
normal derivatives are related by a 2 x 2 transfer matrix: 



(6.20) 



Note, that the transfer conditions (6.20) do not assume identification of cj)'^ and 
(f)~ on E. In other words, there is no ad hoc relation between the restrictions of 
the vector bundles V~^\^ and We can even consider the situation when 

we have dimV^+ ^ dimV"^, i.e. the fields on M~ and M~ can have different 
structures with respect to space-time and internal symmetries. S*^^ are some 
matrix valued functions on E (one can even consider the case when they are 
differential operators). It is clear from the notations on which spaces they act. 
For example, 5*+^ : V^\yi V^^|s- 

The matching conditions (6.20) arise in heat transfer problems [119], some 
problems of quantum mechanics [256], and in conformal field theory [40]. In a 
formal limiting case S++-S-+ = S—-S+~ oo while v = 2{S++ + S+-) is 
kept finite one arrives at transmittal boundary conditions (6.11), (6.14). The 
heat kernel coefficients are divergent in this limit. 

In a particular case of spherical S the heat kernel expansion with transfer 
boundary conditions was evaluated in [81]. General expressions for with 
A; = 0, 1, 2, 3, 4 were obtained in [247]. Somewhat surprisingly, the calculations 
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for (6.20) are easier than for the singular particular case (6.11), (6.14). 

Wc should note that not all singular limiting cases of (6.20) are described 
by the transmittal conditions (6.11), (6.14). The heat kernel expansion for a 
generalisation of transmittal condition is known in the spherically symmetric 
case only [81]. Very little is known about the heat kernel if the transfer matrix 
contains differential operators on E (conformal walls of ref. [40] belong to this 
class of problems). 

The case when the singular potential is located on a surface of co-dimension 
larger than one (i.e. when dimE < n — 1) is rather complicated. Even a care- 
ful translation of this problem to the operator language was done only in 
1960's [63]. Direct calculations show [5,390] that the heat kernel asymptotics 
may contain very unusual (In t)~^ terms. More references can be found in [114]. 

6.4 Non-smooth boundaries 

Rectangular region in a plane is probably the simplest manifold with 
boundaries as far as eigenvalues of the Laplacian are concerned. However, the 
formulae (5.29) - (5.33) are not valid for this case because of the presence 
of corners. The heat kernel expansion on manifolds with piecewise smooth 
boundaries was considered by Kac in his famous paper [286] . He demonstrated 
that for a region in each corner with the inside facing angle a contributes 

, , TT^ — , , 

a2 (corner) = — — (6-21) 

to the coefficient 02 while Oq and Oi are still defined by their "smooth" expres- 
sions. 

The formula (6.21) looks similar to the contribution of a conical singu- 
larity (6.7). The reason for this similarity is that the cone can be obtained 
from the wedge by gluing the sides together and imposing the periodicity 
conditions. 

The study of boundary discontinuities was continued by Apps and Dowker 
who calculated the coefficients 03 [163] and 04 [17] for piecewise smooth bound- 
aries. We also refer to [23] where functional determinants on simplicial com- 
plexes were analysed, and to [161] where divergences in the Casimir energy 
found in [338] were attributed to non-smoothness of the boundaries. A recent 
study [340] should also be mentioned. 

When the angle a goes to we obtain a cusp. In this limit (6.21) is di- 
vergent. Presence of the cusp is an essential singularity which modifies powers 
of the proper time t which appear in the asymptotic expansion of the heat 
kernel [392]. 

6.5 Dielectric bodies 

Calculations of the Casimir energy of a dielectric body have attracted 
much attention and created many controversial results. A (rather large) liter- 
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ature on this subject is reviewed in [79.322]. Quantum field theory formulation 
of this problem is known for a long time already (cf. [380,323, 116]). However, 
the heat kernel analysis of divergences in the Casimir energy in dielectric is a 
relatively new subject. 

Wave propagation with variable speed of hght c{x) is described by the 
operator D = — c(a:)^V^ + . . . , where we neglected the lower order terms. In a 
dielectric medium c is expressed in terms of the dielectric permittivity e and of 
the magnetic permeability /i: c(a;)^ = l/{s{x)iJ,{x)). For a smooth distribution 
of £ and n the operator D is a particular case of curved space Laplacian with an 
effective metric defined by c{x) . The heat kernel coefficients can be calculated 
in the standard way [78]. 

Consider now a dielectric body bounded by S. Typically, c{x) (and the 
effective metric!) jumps on S. This singularity is much stronger than the one 
considered in sec. 6.3. Thus the geometric interpretation of this problem is very 
difficult. Very little is known about the heat kernel expansion in a dielectric 
body of an arbitrary shape. There are calculations for a dielectric ball [77] and 
for a dielectric cylinder [80]. These calculations exhibit a puzzling property of 
the heat kernel expansion in dielectrics in dilute approximation (e ~ 1): for a 
dielectric body the Casimir energy in the ultra violet limit behaves better than 
for the "smooth" case. The heat kernel expansion for a frequency dependent 
e was considered in refs. [202,76]. 
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7 Anomalies 



The most immediate application of the technique developed in the pre- 
vious sections is calculation of quantum anomalies which are defined as (non- 
zero) variations of quantum effective action with respect to symmetry trans- 
formations of the classical theory. In this section we consider two most impor- 
tant examples of quantum anomalies. These are conformal (scale) and chiral 
anomalies. We also discuss briefly the Index Theorem. 

There exists also a broader view on the anomalies which includes any qual- 
itatively new phenomena of a quantum theory which are absent in its' classical 
counterpart. An example is the so-called dimensional reduction anomaly [216]. 
For two operators Di and D2 the quantity det{Di) det{D2) / dct{DiD2) is 
called the multiplicative anomaly [187] since det{DiD2) 7^ det(£'i) det{D2) 
is an "anomalous" property of infinite dimensional operators. It is interesting 
to note that the heat kernel is also useful for treatment of these non-standard 
anomalies. 

To simplify the discussion in this section we work on manifolds without 
boundaries. 



7.1 Conformal anomaly 



Conformal invariance is one of the symmetries which are usually bro- 
ken by quantisation. This phenomenon (called conformal or trace anomaly) is 
known since mid 1970's [118,144,165,96,398] (see [69,174] for more extensive 
literature) . 

The vacuum polarisation induced by quantum effects is described by the 
energy-momentum tensor 

T (71) 

where W is the quantum effective action calculated on the background with 
the metric tensor g^i,. Consider the conformal transformation 

9,u ^ e'P^^^g,, (7.2) 

for an infinitesimal value of the parameter: (yf^j, ^ (1 + 25p)g^v, g^^ ^ (1 — 
25p)g^'' . According to (7.1) variation of the effective action reads: 

5W = \j rx^T.Jg'^'' = - [ rx^T^6p. (7.3) 

It is clear from the equation above that the trace of the energy-momentum ten- 
sor measures conformal non-invariance of the theory. If the classical action C 
is conformally invariant, classical energy-momentum tensor is tracclcss. How- 
ever, even in this case conformal invariance is typically broken by quantum 
effects. For this reason, quantum Tj^ is called trace or conformal anomaly. 
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If the classical action L is conformally invariant, the fluctuation operator 
D is conformally covariant. This means that D transforms homogeneously 
D e-^PD under (7.2) . 

We restrict ourselves to the one-loop level and employ the zeta function 
regularization in which the effective action is expressed through the zeta func- 
tion of the operator D (see (2.32)). Hence, we have to study the conformal 
properties of C,{s^D). One can prove that the variation of the zeta function 
with respect to variation of the operator D reads 

5C{s,D)^-sl}^{{5D)D-^-'). (7.4) 

Next wc use that under infinitesimal conformal transformations the op- 
erator D transforms as 

5D = -2{5p)D . (7.5) 

This equation yields 

5as,D) = 2sas,5p,D). (7.6) 

For the operators which we consider in this section the zeta function is regular 
at s = 0. Consequently, the variation of the effective action (2.32) reads: 

m = -C(0, D) = -a„((5p, D) , (7.7) 

where we have used the relation (2.27) between the zeta function at zero 
argument and the heat kernel coefficient a„. 
We compare (7.3) with (7.7) to see that 

Ti:{x) = an{,x,D). (7.8) 

Note, that C(0, D) (as well as the pole term in (2.31)) is conformally invariant. 
Consequently, the conformal anomaly is not divergent, and (7.8) does not 
contain the normalisation scale ji. 

The anomaly (7.8) is defined by the same coefficient as the divergent 
part of the effective action in the zeta function regularization. Important dif- 
ference is that the divergence is given by the integrated coefficient a„(£)), while 
the conformal anomaly is defined by the localised coefficient an{x, D). The use 
of the integrated coefficient in (7.8) is insufficient to recover total derivatives 
in the anomaly. 



° Strictly speaking this property holds up to a similarity transformation D — > 
e°'PDe^°'P whieh does not change the functional determinant. 

To derive this formula rigorously we have to use the method of [21,368]. We first 
assume that s is sufficiently large to keep us away from the singularities, then use 
the Mcllin transformation (2.24), perform the variation, then perform the transfor- 
mation back, and then continue the result to s = 0. This is a perfectly standard 
procedure which allows us to work with variations of positive integer powers of D; 
see [198] for further detail. 
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As an example, consider qTiantum scalar field in two dimensions'^ 
coupled to the background dilaton 

jC^J d^x^e-^''{d,^){dA)9'^ . (7.9) 

On dimensional and symmetry grounds the inner product may also contain 
an arbitrary function ^ of the dilaton: 



J d^x^e-^''Mx)Mx)- (7.10) 



Such couplings and inner products appear, for example, after the spherical 
reduction of higher dimensional theories to two dimension (see the review 
paper [265] for more details and further references). 

The rescaled field (p = e~^(j) possesses the standard dilaton-independent 
inner product (3.3). In terms of the action (7.9) reads 

£ = y (fx^4>D4>, (7.11) 

D = -e'(*-*)9^^(V^V. + 2{^>,, - $,^)9, 

+*^, + ^,^^,,-2^,^$,), (7.12) 

where comma denotes covariant differentiation with V: = Vj^^. We can 
bring D to the standard form (2.2) 

D^-{gi^-V^V, + E) (7.13) 

by introducing the effective metric g^'^ = e^^^~'^^g^'^ and the covariant deriva- 
tive 

V^ = d^ + t^ + uj^, o;^ = - , (7.14) 

where V is the Christoffel connection for the metric g. Here the potential E 
reads 

^ = ^^'^(-^,A^ + ^,M.)- (7.15) 
Now we combine (7.8) with (4.27) and the definitions given above to obtain: 

r/;' = ^(i?-6(V*)2 + 4V2* + 2V2*). (7.16) 

In the case \1/ = $ the expression (7.16) was first obtained in [336], and in the 
general case in [303] . 



Much work on conformal anomalies on two-dimensional manifolds was done in 
the context of string theory. Here we like to mention the papers [10,215,214,113,11]. 
'^Different values for the numerical coefficients in the conformal anomaly (7.16) 
were reported in [86,344,321,345]. The reason for these discrepancies was clearly 
stated by Dowker [160] who confirmed the result (7.16). More extensive literature 
can be found in [265]. 
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In four dimensions the conformal anomaly for different spins can be read 
off from (4.35) with the numerical coefficients given in Table 1. 



1.2 Chiral anomaly 

Chiral anomaly was discovered in 1969 by Adler, Bell and Jackiw [3, 60] 

and since that time plays a crucial role in understanding of the low energy 
hadron physics. A detailed introduction to the field and extensive literature 
can be found in [65,428]. 

The spinor action (3.22) with the Dirac operator given in (3.26) is invari- 
ant under the gauge transformations 

6xij = -X^j , dxip = ipX , 
6,Al = [Al A] , 

6xA^ = d^X+[A^,X] (7.17) 
and local (Euclidean) chiral transformations 

~S^Al = d^if + [A^, if] , 

6^A, = -[Al^] (7.18) 

with anti-hermitian local matrix parameters A and (p. The Dirac operator 
transforms as: 

Sxl/)^\]p,X], ~Sj)^tWj',p}. (7.19) 

The Lie algebra structure of the transformations (7.17) and (7.18) is encoded 
in the following relations: 

hih2 - ^2^1 = 5[Ai,A2] (7.20) 
- S^ph = S[x,ip] (7-21) 

^ipi^ip2 ~ ^ip2^^Pl ~ ~^lipi,f2] ) (7.22) 

where all transformation parameters are taken at the same space-time point. If 
the matrices A(x) and ^p{x) belong to a finite-dimensional compact Lie algebra 
of some Lie group Q, the transformations (7.17) and (7.18) generate locally 
the group QxQ with gauge transformations belonging the diagonal sub-group. 

The gauge invariance can (and should) be retained in quantum theory 
while the chiral invariance is typically broken by the quantisation. Let us 
study these effects in the zeta function rcgularization. Define determinant 
of the Dirac operator as a square root of the determinant of the associated 
Laplace operator: 

log det^ = ^ log det|>^ = ^ log det D (7.23) 
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The effective action W corresponding to the path integral (3.41) reads 

W^-logZ^-^ logdet D = ^C(0, D)' + ^ ln(/x2)C(0, D) , (7.24) 

where // is a normahsation scale. 
By virtue of (7.4) 

S^Cis, D) = - (sTr(p, \]D-^-')) = (Tr(p-^ A])) = . (7.25) 

This proves gauge invariance of the effective action. 
For the chiral transformation we have: 

A{ip) := S^W = -2iTr(7V/^-1|s=o • (7.26) 

A{ip) measures non-invariance of the effective action with respect to the chiral 
transformations. It is called the chiral anomaly. One can express A through 
the heat kernel coefficients: 

A{ip)^-2tan{-f'ip,D), (7.27) 

where an{'y^(p,D) is defined as in (2.21) but with a matrix-valued smearing 
function instead of the scalar one. 

The chiral anomaly A{f) should satisfy certain consistency conditions 
following from the Lie algebra identities (7.20)-(7.22) and gauge invariance of 
the effective action: 

SxA{^)=A{[\,^]), (7.28) 
~S<p,AM-K2A'Pi)-Q- (7-29) 

The relations (7.28) and (7.28) are called the Wess-Zumino consistency con- 
ditions [417]. One can check by a direct calculation that the anomaly defined 
in this section indeed satisfies these conditions. For this reason A is called the 
consistent anomaly (as opposed to the covariant anomaly which we do not 
consider here). 

Let us turn now to calculation of the chiral anomaly A. Since the smearing 
function 7^</7 in the heat kernel coefficient (7.27) is matrix- valued we need more 

information than we possess at the moment. To recover the missing terms 
one can adopt the strategy of the paper [94]. Consider Tt{Q exp(—tD)) with 
arbitrary matrix-valued function Q and arbitrary Laplace type operator D. 
There is an asymptotic expansion 

TrL2{QeM-tD)) ^ Y.^^'^'^^'^'kiQ, D) (7.30) 

fc>0 

where the coefficients Cfc are locally computable. This means that they can 
be represented as integrals of local invariants constructed from Q and local 
invariants of the operator D. These local invariants enter with some coefficients 
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(universal constants) which are to be defined. Since Q docs not commute with 
E and Q there are more invariants than before. In the particular case Q = lyf 
(where ly is the unit matrix and / is a function) we should recover the old 
result (4.26) - (4.29). For k — 0,2,4 this last requirement is strong enough 
to recover ak{Q,D) completely. Therefore, calculation of the chiral anomaly 
in the dimensions n = 2 and n = 4 is quite simple: take the expression 
(4.27) or (4.28) and substitute 2^7 V for / and (3.27), (3.28) for u^,, E and 
Vt^i,. Computation of some 7-matrix traces is still required, but nevertheless 
calculations are considerably easier than presented in many papers. In flat 
space and with — Q the result is particularly simple: 

n = 2 : A{^) = / d^xtr{ipe'"'F^,) (7.31) 

n = 4 : A{ip) = Y^I d^xtr{ipe'''"'''F^,F,,) (7.32) 

where tr denotes trace over internal (flavour or colour) indices. We can easily 
generalise this resTilt for arbitrary even dimension n. First we observe that the 
only way to construct a pseudoscalar density of appropriate dimension (for 
= 0) is to contract n/2 tensors F^^, with the Levi-Civita tensor e^^^^---^"^. 
Such tensor structure can be produced only by a trace of 7^ and the maximal 
number {n) of the gamma matrices. The only invariant having the required 
form is 

j (T xtriQE''/^) . (7.33) 

This term does not contain derivatives. Under the trace Q and E commute. 
Therefore, we can calculate the coefficient in front of (7.33) by considering the 
particular case £' ~ x const.. For this simple case the dependence of the 
heat kernel on E is given simply by e*^. By picking up an appropriate term in 
the expansion of the exponential we find that the coefficient in front of (7.33) 
in a„ is (47r)""/^((n/2)!)"\ Next we substitute E = \^'^YF^,v and take trace 
over the spinorial indices to obtain: 

^^^^ = 1 ^""^'^ (^^""■■■"^-- • • • ) ^^-34) 

For global chiral transformations {ip = const.) this formula was obtained in 
[421]. 

As soon as we know the anomaly (7.34) for = 0, many more terms 
(containing A^) can be restored by using the Wess-Zumino condition (7.29). 

Some comments on the Fujikawa approach [220,221] to the chiral anomaly 
are in order. Consider the path integral (3.41). Since the action is invariant 
under the chiral rotations, the only source of the chiral anomaly could be non- 
invariance of the path integral measure. The Jacobian factor appearing due 
to the change of the spinorial variables (7.18) can be formally represented as 

J = det(l - 2i(p-f^) . (7.35) 
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Then, to the first order of </?, 

A{ip) = 6^W = - log J-^ ~ -2iTr(937^) (7.36) 

where has appeared due to the negative homogeneity of the fcrmionic 
measure. The operator Lp^^ is not trace class on the space of square intcgrable 
spinors (though formally its' trace is zero at every point). The right hand side 
of (7.36) is therefore ill defined. Fujikawa suggested to replace (7.36) by a 
regularized expression 

A(ip) = -2i lim Tr (cp-f^e-^/^") . (7.37) 

M— >oo ^ ' 

The chiral anomaly (7.37) has now the heat-kernel form with the identification 
t — M~^. If we suppose that all positive powers of the regularization parameter 
M in the small t (large M) asymptotic expansion are somehow absorbed in 
the renormalization, we arrive at 

A{ip) = -2m„(7V, D) (7.38) 

that is just the expression (7.27) obtained above in the zeta function regular- 
ization. This method, based on calculations of the regularized Jacobians, can 
be apphed to conformal anomaly as well [222] . 

The most essential ingredient of the anomaly calculation presented in this 
section is the homogeneous transformation low (7.19) for the Dirac operator. 
This homogeneity allowed us to restore the power —s in the transformation law 
(7.4) for the zeta function and to obtain a simple expression for the anomaly 
(7.26). It is clear therefore that as long as the operators transform homoge- 
neously we shall obtain relatively simple local expressions for corresponding 
anomalies (understood as variations of the effective action with respect to 
infinitesimal transformation of the background fields). This suggests to con- 
sider extensions of the chiral group. For example, one of such extensions [346] 
identifies the group parameters with the diquark fields. 

7.3 Remarks on the Index Theorem 

This report is mainly devoted to local aspects of the heat kernel expansion. 
There is however one "global" apphcation of the heat kernel which cannot be 
ignored. This is the Index Theorem. In this section we briefly sketch formal 
mathematical aspects of the index construction and its' relation to the heat 
kernel. For more details we refer to [177,245,189,65,428]. Physical applications 
of the index theorem to gravity, gauge theories and strings are so numerous 
that we cannot even mention them all. 

The index theorem was first formulated by Atiyah and Singer [22], and 
the heat kernel approach appeared later in [236,18]. Roughly speaking, their 

Applications of this method to topological anomalies can be found in [223] . 
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construction is as follows. Consider two vector bundles Vi and V2 over a man- 
ifold M. Let the operator P map Vi to V2, and let \4 2 have non-degenerate 
inner products ( , )i 2 so that one can define an adjoint P"^ by the equation 
(02,-P0i)2 = {P^4>2,<Pi)i- As an example, one can keep in mind Vi (respec- 
tively, V2) describing positive (respectively, negative) chirality spinors. In this 
example P is a part of the Dirac operator (see eq. (7.44) below). Let us define 
two operators Di = P'^P and D2 = PP'^ acting on (smooth sections of) Vi 
and V2, and let us suppose that Di and D2 are elliptic. In this case we deal 
with an elliptic complex. 

Since elliptic operators may have only a finite number of zero modes 
(say, A^i and A^2 for Di and D2 respectively), we may define the index by the 
following equations: 

index(P) = Ni-N2 

— dim ker Di — dim ker D2 

= dim ker P - dim ker P^ . (7.39) 

We also have the intertwining relations 

P>iPt = ptp)2 , PP*! = P>2P , (7.40) 

which tell us that non-zero eigenvalues of Di and D2 coincide. Consider now 
the heat kernels, 

K{t, Di) - K{t, D2) = e"*^' - e"*^' = A^i - ^"2 = index(P) (7.41) 

(Ai) (A2) 

(where Ai_2 are eigenvalues of 1^1,2) ■ 

To be more specific, let us suppose that Di and D2 are of Laplace type and 
that boundary conditions and singularities (if any) are such that the power- 
law expansion (2.21) exists. Then, expanding both sides of (7.41) in a power 
series of t one obtains: 

ak{Di) - ak{D2) = Q for k ^ n, 

an{Di) - an{D2) = index(P) . (7.42) 

The last equation provides a simple way to calculate the index from known 
heat kernel expansion. It also allows to understand a very important property 
of the index, namely, that it is a homotopy invariant under quite general as- 
sumptions. Indeed, suppose that P depends on a parameter a in such a way 
that all geometric quantities (including the metric, connection, the matrix po- 
tential, the boundary conditions etc.) corresponding to Di and P>2 are smooth 
functions of a. Clearly, no "essential" deformations, like changing the order of 
the operator, or adding higher derivative terms to the boundary conditions, or 
including new types of the singularities, are allowed. Under these smoothness 
assumptions a„(Di) and a„(D2) are smooth functions of a as well. Therefore, 



69 



the index is also a smooth function of a. Since the index is an integer, it 
can only be a constant. Hence, index(P) is invariant under the deformations 
described above. 

The relations (7.42) are also useful for calculations of the heat kernel 
coefficients as they give restrictions on the universal constants appearing in 
front of independent invariants (see, e.g., [245]). 

Let us consider an example of the spin complex. Let M be an even di- 
mensional manifold admitting a spin structure (so that one can define spinor 
fields on M) . In a suitable basis the chirahty matrix may be presented in the 
diagonal form (at this point it may be useful to consult sec. 3.3): 



75 = . (7.43) 




Let Vi correspond the +1 eigenvalue of 75 (positive chirality), and V2 ^ to 
— 1. The Dirac operator (3.26) anti-commutes with 75 and, therefore, may be 
represented in the form: 

^ pt' 



P 



" I • (7.44) 

We also have: 

# = 1 _ 1 • (7-45) 




We conclude that index(P), which now measures difference between dimen- 
sions of the null subspaces of p with positive and negative chirality, is equal 
to global chiral anomaly: 

index(P) = IV (e-*^^^) - IV (e"*^^^) = IV (756-*^') = a„(75,^') (7.46) 

Therefore, the index can be calculated by the methods of the previous section. 
In particular, the index of the Dirac operator in two and four dimensions on 
fiat manifolds without background axial vector fields (A^ = 0) can be read off 
from (7.31) and (7.32) with cp :— i/2. We also see that chiral anomaly provides 
very important information on topology of the background. 

One can associate an index with any elliptic complex which must not 
necessarily contain only two vector bundles. In general, one has a sequence 
of bundles Vp and a family of operators Pp which map (smooth sections of) 
Vp to (smooth sections of) Vp+i. One requires that PpPp+i = and that 
the corresponding "Laplacians" Dp = P^Pp + Pp_iP^_i are elhptic. Important 
examples include the dc Rham complex (for which Vp consists of p-forms on M 
and P is the exterior derivative), the Dolbeaux complex which deals with the 
forms on a complex manifold, and the signature complex which treats selfdual 
and anti-selfdual forms. More details can be found in [177]. Another very 
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important construction of this type is the Witten index [422] which constraints 
the supersymmetry breaking. 
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8 Resummation of the heat kernel expansion 

The heat kernel coefficients define the one-loop counterterms in the back- 
ground field formalism. In many cases the heat kernel can also give a useful 
information on the finite part of the effective action. Just one of the examples 
is the large mass expansion (1.22) which is valid when all background fields 
and their derivatives are small compared to the mass of the quantum field. In 
order to get the effective action in other limiting cases one has to re-arrange 
the heat kernel expansion . 

8.1 Modified large mass expansion 

In many physical applications there is a quantity Ai^ which is large com- 
pared to the rest of the background fields and their derivatives. Therefore, it 
is a well motivated problem to construct an expansion of the effective action 
in a power series of M.~^. To do so, one has to re-express the heat kernel as 

K{t; D) = e-'^' ^ t^''-''^'^~ak . (8.1) 

k 

To obtain the effective action one has to integrate the heat kernel over t (cf 
eq. (1.18)). To simplify the argumentation we assume the cut-off regulariza- 
tion (1.20), although in other regularization schemes the result will be essen- 
tially the same. The divergent and finite parts of the effective action are given 
by (1.21) and (1.22) respectively up to the obvious replacements: m Ai, 
{ATi)~^f%k — * Sfc. After a suitable renormalization, this procedure indeed gives 
a large M. expansion of the effective action. 

Therefore, the problem is reduced to calculation of a^. If Al commutes 
with D the coefficient are the heat kernel coefficients for the operator 
D — Ai^. This case returns us to the standard large mass expansion. If Ai^ 
does not commute with D then 

To achieve an equality the right hand side of (8.2) must be corrected by 
commutator terms. In this case calculation of requires some amounts of 
extra work. 

Heavy particles of non-equal masses described by the mass matrix A4^ 
is probably the most immediate example of a physical system to which the 
modified large mass expansion should be applied. Corresponding technical 
tools were developed recently [355,354,353]. 

The next example is a scalar field in curved space (cf. sec 3.1). Parker 
and Toms suggested [359, 358] to use the modified large mass expansion with 

A diagrammatic technique which can be used in resummations of the heat kernel 
expansion is described in [331]. 
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A4'^ = + which partially sums up contributions of the scalar curvature 
R to the effective action. This formalism was developed further in [284,393]. 

In some cases the term-by-term integration of the heat kernel expansion 
gives good estimates of the vacuum energy even if no large parameter is ex- 
plicitly present in the model (see, e.g. calculations of quantum corrections to 
the mass of two dimensional solitons [176,9]). 



8.2 Covariant perturbation theory 



Suppose that the matrix potential E, the Riemann curvature, and the 
field strength Q, are small but rapidly varying. Can we get any information 
on the heat kernel coefficients containing a fixed power of the quantities listed 
above and arbitrary number of derivatives? For the linear and quadratic orders 
on a manifold without boundaries the answer is positive and the results can 
be obtained either by the functorial methods of sec. 4.1 [239,92] or by solving 
the DeWitt equation [24,25,26] (cf. sec. 4.3). Retaining only the leading terms 
we have for A; > 1: 

a2k{x, D) = (47r)-"/2 \ai{k)M'-^R + a2{k)M'-^E (8.3) 



where [237, 239] 



-l-higher order terms] , 



One can sum up the expansion (8.3) and corresponding higher order terms 
to obtain an information on the behaviour of the full heat kernel in the limit 
described above. There is, however, a more straightforward method to control 
non-localities which is called the covariant perturbation theory [56,57,58]. To 
make the idea of this method most transparent we consider a simple example 
of flat M and zero connection u. The exponent exp(— iD) = exp(i(A -|- E)) 
with A = 9^ can be expanded in a power series in E: 



S2 



J dse^'-'^^Ee'^+ J ds2j dsie^'-''^^Ee^''-''^^Ee''^+. . . (8.5) 



The heat trace can be also expanded, 

oo 

X(t,D) = Tr(e-*^)=E^.W, (8.6) 

j=0 

where Kj contains the j'th power of E. 

Covariant perturbation theory approach prescribes to take the 0th order 
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heat kernel in the free space form (cf. eq. (1.12) for m — 0): 

Ko(x, y; t) = (Ant)-"/' exp (-^^^) , (8-7) 

which is an exact kernel on M = M"^ only^^ . This formula neglects all global 
contributions and, therefore, is valid only for sufficiently close x and y and for 
small t. We have: 

Ko{t) = / d'^x^trvKoix^x-t) = (Ant)-''/' [ d^'x^trvilv) , (8.8) 

JM JM 

where trylv simply counts discrete (spin and internal) indices of D. This 
formula reproduces the ao contribution to the heat kernel. 
In the next order, we have: 



K,{t) = Tr dse^'-'^'^Ee'^^ = Tr (^J^ dse'^E^^ 
= nry / d"'x^/gKo{x,x■,t)E{x) 

JM 

where we used the cyclic property of the trace and the expression (8.7) for 
the 0th order heat kernel. This expression is consistent with the k — 1 term 
of (8.3). The terms with k > 1 are total derivatives, and, therefore, they do 
not contribute to the integrated heat kernel. 

The quadratic order of the heat kernel reads 

.52 



K2{t) ^TvUdsij dsie^'-''^^Ee^''-''^^Ee''^ 



\G 

t S2 



= try / dy f dz ds2 dsiKo{z,y;t-S2 + Si)E{y) 

JM JM J J 


xKo{y,z;s2-si)E{z). (8.10) 

Next we introduce the rescaled variables i = s/t and get rid of redundant 
integrations to obtain: 

2 ^ 

K2{t) = %iivl dy [ dz I diKo{z,y-t{l-i))E{z) 

I JM JM J 


~KK^{y,z-ti)E{z). (8.11) 



The original paper [56] contained a curved space generalisation of (8.7). This, 
however, does not improve the global issues discussed below. 
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Now we use the identity 

Ko{z, y- t(l - 0)Ko{y, t^) = i4nt)-'^/^Ko{z, y; t^l - 0) (8-12) 

and relate the heat kernel on the right hand side of (8.12) to a matrix element 
of exp(t^(l — ^)A) to obtain the final result: 

m) = Jj^^x^Ef{-tA)E, (8.13) 

where the non-local form-factor / reads: 

m = \] d^e-'^^^'-^^ = ie-^/^^^AErfi [^/2] . (8.14) 



As we have already discussed above, applicability of this formula is limited by 
our choice of KQ{x,y;t). Namely, the potential E{x) should have a compact 
support and t should be reasonably small. Note, that the small t approximation 
in the expansion (8.5) is self- consistent: if t is small, the integration variables 
Si are even smaller. 

The main difficulty in constructing an expansion in powers the Riemann 
curvature and of the field strength is to organise the procedure in a covari- 
ant way. The details of the construction and higher order form-factors can 
be found in [56,57,58,50,49]. From further developments of this method we 
mention the work of Gusev and Zelnikov [267] who demonstrated that in two 
dimension one can achieve considerable simplifications in the perturbation ex- 
pansion by using the dilaton parametrisation of the potential (cf. eq. (7.15)). 
Recently Barvinsky and Mukhanov [53] suggested a new method for calcula- 
tion of the non-local part of the effective action based on the resummation of 
the perturbation series for the heat kernel. This method was extended in [48] 
to include late time asymptotics of the heat kernel in curved space. 

Let us stress that the expansion (8.5) can be used also for singular poten- 
tials. For example, it is very effective for the calculation of the heat kernel for 
5-potentials concentrated on a co-dimension one subsurface [231,82,330] (cf. 
sec. 6.3). If the 5-potential has its support on a submanifold of co-dimension 
greater than one, the expansion diverges [82]. With a suitable choice of the 
zeroth order heat kernel Ko{x,y;t) and of an operator to replace E in (8.5) 
one can treat manifolds with boundaries [74] where the perturbative expansion 
takes the form of the multiple refiection expansion [42,43,272]. 

8.3 "Low energy" expansion 

Let us now turn to the opposite case when the derivatives are less impor- 
tant when the potential and the curvatures. In this case one has to collect the 
terms which are of a fixed order in the derivatives, but contain all powers of 
E, R and Jl. Since the derivatives are sometimes identified with the energy. 
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this approximation is being called the low energy expansion. This scheme goes 
back to Schwinger's calculations [379] in constant electromagnetic fields . 

Let us consider a simple example [71]. Let M = M^. Consider a scalar 
particle in constant electromagnetic field with the field strength F12 — — -F21 = 
B. As a potential we choose Ai = 0, = Bx^. Then the operator acting on 
quantum fiuctuations is 

D = -dl-{d2-iBx^ f + . (8.15) 

This operator commutes with 82- Therefore, we can look for the eigenfunctions 
of D in the form = exp(iA;x^)0(x^). 

D(i)k{x) = (^-dl + B\x^ - k/Bf + rn^) ^^(a;) . (8.16) 

In the x^ direction we have the one- dimensional harmonic oscillator potential 
(cf. (6.3)). Therefore, the eigenvalues are A^^p = {2p + 1)|-B| + m^. These 
eigenvalues do not depend on k. For this reason the heat kernel 

X(i,L>)=^exp(-afe,p) (8.17) 

k,p 

is ill defined. To overcome this difficulty it was suggested [71] to put the system 
in a box (without specifying any boundary conditions, however) and replace 
the sum over k by the degeneracy factor (|i?|Vol)/(27r), where Vol is volume 
of the box. The degeneracy factor is chosen in such a way that the resulting 
heat kernel 

Kit, D) = ^^e-™'*[sinh(5t)]-^ (8.18) 

4:71 

behaves as Vol/(47rt) in the limit i? — > 0. 

Calculations of the effective action in covariantly constant background 
gauge fields by the spectral theory methods have been performed by many 

authors [59,151,185,183,71]. Such calculations were motivated, at least par- 
tially, by various models of quark confinement (see [356, 4] and references 
therein) . 

There exists an algebraic method [27,29] which allows to evaluate the low 
energy heat kernel by using exclusively the commutator algebra. Let us briefiy 
formulate the results of [27, 29] . Consider a flat manifold {Rnvpc — 0) with the 
background fields satisfying the "low-energy conditions" : 

V^Q^, = 0, V^V,Vp£; = (8.19) 

with usual definitions of E and Vt (see (2.1) - (2.10)). Moreover, let us suppose 
that the background is "approximately abelian", i.e. that fi^^,, E and all their 
covariant derivatives are mutually commuting matrices. 

We like to mention also the paper [275] which treated the effective action in 
external electromagnetic field from a different point view. 
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With all these assumptions a closed expression for the heat kernel may 
be obtained: 

K{t; X, x; D) = {AT^ty'^ exp i^E + + jE.^^'"' E.}j , (8.20) 

where $ and ^ are comphcated functions of t, E and [29]. If V^V,^-E = 0, 
^ , , 1 , , / sinhin\ 

^{t) = {my^ {mcoth{m) - 1) , (8.21) 

where Q has to be understood as a space-time matrix Q^^, so that multipli- 
cation in (8.21) is the matrix multiplication, and det is the determinant of an 
n X n matrix. These formulae generalise the equation (8.18) and justify the 
choice of the degeneracy factor made to derive it. 
If Q = 0, 

1, 1 /sinh(2t^\ 
$ t = --In det ^ J- , 

^{t)^- (tVPy^ (tanh(iVP) - tVP)) , (8.22) 

where P^j, = — (l/2)V^Vjy-E. In the particular case of one- dimensional har- 
monic oscillator these formulae reproduce (6.5). 

In curved space the best one can do in the framework of the low-energy 
expansion is to consider locally symmetric manifolds, i.e. vanishing fi, co- 
variantly constant Riemann curvature and constant E. In this case formulae 
similar to the ones presented above are available [28,30]. Covariantly constant 
curvature means that locally the manifold M is a symmetric space. Various 
approaches to the heat kernel on such manifolds are described in detail in 
monographs and survey articles [282,115,109]. In particular, very detailed in- 
formation may be obtained for group manifolds (see, for example, [153]) and 
for hyperbolic spaces [110]. 

8.4 Heat kernel on homogeneous spaces 

In this section we briefly explain how one can find the spectrum of some 
"natural" differential operators on homogeneous spaces by purely algebraic 

methods. We start with some basic facts from differential geometry and har- 
monic analysis [301,276]. Consider a homogeneous space G/H of two compact 
finite-dimensional Lie groups G and H. The Lie algebra ^ of G can be decom- 
posed as 

g^n®M, (8.23) 
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where V. is the Lie algebra of H and M. is the complement oiTi in Q with 
respect to some bi-invariant metric. We have: 

[n,M\^M, (8.24) 

where [ , ] is the Lie bracket on Q. If, moreover, [A^, A^] C Ti, then G/H is a 
symmetric space. We do not impose this restriction. 

M. can be identified with the tangent space to M — G/H ai the origin 
(i.e. at the point which represents the unit element of G). Eq. (8.24) tells us 
that H acts on M. by some (orthogonal) representation. This action defines 
the embedding 

n C so{n) . (8.25) 

All physical fields are classified according to certain representations of the Lie 
algebra so{n). Restrictions of these representations to H define transformation 
properties of the field with respect to Tl. Prom now on we work with each 
irreducible representation of Ti separately. 

The field belonging to an irreducible representation T[H) can be 
expanded as (see, e.g., [373]): 

^a{x) = (Vol)-i I¥'^AL{9-.')<t>^l , (8.26) 

where Vol is the volume of G/H, dr = dimT(iJ). We sum over the repre- 
sentations T^^^ of G which give T[H) after reduction to H. ^ labels multiple 
components T[H) in the branching T^^ [ H, dj — dimT^.'^ q runs from 
1 through dj. The matrix elements of T^^^ have the following orthogonality 
property 

/ d''x^T^ll{g;')T^^l{g;') = (Vol)rf7^rfr%5p/,y (8.27) 

Therefore, to construct the harmonic expansion on G/if it is necessary to have 
powerful methods for reduction of the representations from G to H. There are 
several standard [46,387] and less standard [311,310] techniques which may 
be used depending on the particular homogeneous space. 

It is important that not only the harmonic expansion but also the spec- 
trum of the invariant operators can be analysed by the group theoretical meth- 
ods. The covariant derivative on G/H reads [301]: 

= + rjf] . (8.28) 

Here V''^' is the canonical covariant derivative on G/H. At the origin 
can be identified with the tangent space generators from Ai taken in the 
representation T^^\ The part r|f^l depends on the invariant metric on G/H 
and on the structure constants of ^ restricted to M.. On symmetric spaces such 
structure constants are zero and, therefore, the Laplace-Beltrami operator has 
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a particularly simple form: 

D ~ V^t^lvJ,^^ ~ C2{G) - C2{H) , (8.29) 

where C2 are quadratic Casimir operators of G and H which depend on the 
representations T*^*) and T{H) respectively. On general homogeneous spaces 
the expressions are a bit more complicated (see, e.g., [407,310] for explicit 
examples). In any case, eigenvalues of D are given by a second order poly- 
nomial Q{mi, . . . ,mfc) of several natural numbers mi. These eigenvalues are 
in general degenerate with multiplicities defined essentially by dimensions of 
the representations of G and H. They are also polynomials^^ in mi. The heat 
kernel is then represented as an infinite sum 

K{t, D) = Y^ N{mi, . . . , mfc) exp(-tQ(mi, . . . , m^)) . (8.30) 

mi 

There exist several tricks which can be used to evaluate the t — > asymptotics 
of such sums. For example, one may use the Poisson summation formula (4.4). 
By taking derivatives with respect to t one obtains 

^ P |^_^^2^ ^ 0F2-'-(2r - l)!!r(2'-+i)/2 + O (e-^/*) (8.31) 
for r G N. 

More general polynomials of I may be treated by the Euler-Maclaurin 
formula which reduces sums to the integrals. Let F{t) be a function defined 
on < r < 00. If the 2m-th derivative F(^"*)(t) is absolutely integrable on 
(0,00) 

^ roo 1 

Y.F{i)- F{r)dr^-[F{0) + F{k)] 
^ io 2 

m—l TD 

+ E 77^,[F''''-'\k) - F^'^-'\Q)] + Rem„(n), (8.32) 

s=l 

where the reminder Rem^ satisfies: 




(8.33) 



Bs are the Bernoulli numbers. If F{i) is taken to be Ni exp(— tAj), one can take 
limit A; — > 00 in (8.32) and restrict the summation to some finite m to obtain 
asymptotic series for the heat kernel. For example, one can easily recover the 
expansion (8.31). 

The spectrum of differential operators on coset spaces can be calculated exactly 
even if one adds a homogeneous gauge field (see [152] and references therein). 
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In many cases calculation can be done by means of the Mellin transform. 
This method will be clear from the following example. Consider 

oo 

K{t) = ^ A;exp(-tA;^) (8.34) 

k=l 

The Mellin transform of K{t) is 




exp(-xA;2) = ^ A;1-2T(s) = Cr{2s - 1)T (s) , (8.35) 

k=l 



where is the Riemann zeta-function. Inverse Mellin transformation gives 

Kit) = ^ / t-%R{2s - l)ris)ds. (8.36) 

Zm Jc 

The contour C covers all poles of the integrand at s = 1, 0, —1, —2, .... By 
calculating the residues we obtain the desired expansion: 

^^■W 4 - 150 -li-iSo + 

Actually, the techniques introduced above should only be used for com- 
plicated multi-parameter sums (see, e.g., [310]). For a simple one-parameter 
sum it is easier to use combinations of known asymptotic series (cf. Appendix 
of Ref. [67]). 

In this method complexity of calculations of the heat kernel coefficients 
is almost independent of k. Therefore, the algebraic techniques were applied to 
higher dimensional theories where one needs higher heat kernel coefficients to 
perform renormalization or to calculate the anomalies. Rather naturally, the 
most simple toroidal spaces were considered ffist [15, 16,372,302,210, 102], and 
then computations on spheres were performed [117,295,290,67,68,249,105, 
106,180]. Other homogeneous spaces were considered for example in [310,408]. 
In the same approach non-minimal operators on homogeneous spaces were 
treated in [401,403]. We refer to [115] for a more extensive literature survey. 
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9 Exact results for the effective action 

In the previous section we have seen that the heat kernel can be calculated 
exactly for all values of t if the background satisfies certain symmetry or 
smoothness conditions. In this section it will be demonstrated that for some 
classes of the operators the effective action can be calculated exactly without 
imposing any restrictions on behaviour of the background fields. These are 
the cases when the variation of C'(0) with respect to the background can be 
reduced to the heat kernel coefficients which are locally computable. Then the 
variation is integrated to the effective action. 

The most immediate example of the variations admissible for this scheme 
is quantum anomalies (cf. sec. 7). Historically, the first action obtained by 
integration of an anomaly was the chiral Wess-Zumino action. We do not 
consider this action here since such a consideration would require a lot of 
additional technical devices. The reader can consult excellent original papers 
[417, 424] 34 . The next subsect ion is devoted to another well know example, 
which is the Polyakov action obtained through integration of the conformal 
anomaly in two dimensions. 

In sec. 9.2 below we consider a more complicated situation when relevant 
variation of the effective action leads to a linear combination of the heat kernel 
coefficients of several different operators. Typically in this case one obtains 
exact relations between the effective actions rather than the effective actions 
themselves. 

9.1 The Polyakov action 

Let us consider a two-dimensional Riemannian manifold M without bound- 
ary and a scalar field minimally coupled to the geometry. This means that 
we choose $ = in the action (7.9) and \E' = in the inner product (7.10). In 
two dimensions any metric is conformally fiat, i.e. by a suitable choice of the 
coordinates one may transform it to the conformal gauge 

Qi.. = e'^'^^^ri^,, 77^, = diag(l, 1). (9.1) 

In this gauge 

^R=-2r]^'''d,d,p. (9.2) 

The effective action W depends on p only (since there is no other external 
field in this problem). We substitute (7.16) with $ = ^' = in (7.3) to obtain 

iZTT JM 



Calculation of the effective action in two dimensional QED by integration of the 
anomaly may be found in [279]. 
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This equation can be integrated with the initial condition W{p — 0) —0. 

W^^f d^xpri^'^d^d^p. (9.4) 

By returning to the covariant notations we obtain the famous Polyakov action 
[360] (see also [309]): 

Note, that on a compact manifold one has to take into account contribu- 
tions from the zero modes which lead to an additional term in the Polyakov 
action [157]. 

The problem of finding an action whose conformal variation reproduces 
the conformal anomaly can be posed in higher dimension as well [144] and 
some physical consequences of such actions can be studied [218,219]. Com- 
plete expressions for the conformal anomaly induced effective action in four 
dimensions were obtained in [366,100]. Conformal action in the presence of 
boundaries was constructed in [170] . One should keep in mind that unless the 
background is conformally flat the conformal action represents only a part of 
the full effective action. Applicability and limitations of the conformal action 
approach are discussed in [188,41,267]. 



9.2 Duality symmetry of the effective action 



In this section we deal with variations of the effective action which can be 
expressed through zeta functions of several operators (in contrast to genuine 
anomalous transformations of sec. 7 which involve a single zeta function each). 

A rather simple example can be found in the paper [304] which considered 
two operators of Dirac type: 



(9.6) 



on flat two-dimensional manifold without boundary depending on two scalar 
functions $ and \t'. Let us now consider small variations of $ and ^ . According 
to (7.4) the zeta function of corresponding Laplacian changes as 



SC{s,pp^) = -2sTr [(^^^)-M* - {Jp^pyS^ 
Variation of C'(0) reads: 

(5C'(0,#t) = _2 [^(0, Sm,pp^) - C(0, S<^,p^p)] 



(9.7) 



:-2 



(9.8) 



Here we used that pp'^ andp^p are operators of Laplace type, and, therefore. 
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the corresponding zeta functions are regular at s = 0. In the last line we 
replaced ({0) by the heat kernel coefficients by means of (2.27). 

Next we have to rewrite pp^ and p'^p in the canonical form (2.2) with 
the help of (2.3), (2.4). For pp^ we have: 

= i?=-2($-*)/, (9.9) 

where comma denotes ordinary partial derivatives, all contractions are per- 
formed with the effective metric g^'^. Note, that in two dimensions F'^^g^^'^ — 
in the conformal coordinates. To obtain relevant geometric quantities for p^p 
one has to replace $ by — ^ and vice versa. Equation (4.27) gives: 

SC{0,pp^) = I d^x [5^(A^ + 2A$) + 5$(A$ + 2A^)] , (9.10) 

where A is the flat space Laplacian. Let us recall, that C'(0, D) defines the 
functional determinant of D (see eq. (2.33)). The variation (9.10) can be in- 
tegrated to give 

lndet(^^t) = -C'(0,^^^) = — I rf2a;(^A^ + 4^A$ + $A$) . (9.11) 

127r J 

We see from (9.11) that lndet(^^^) = lndet(^^^) which is obviously true 
up to possible topological contributions from zero modes which have been 
neglected in this calculation. The result (9.11) was confirmed in [267] by the 
methods of covariant perturbation theory. 

This example is clearly of (at least an academic) interest since it adds 
up to very few cases discussed above in this section when a closed analytic 
expression for the determinant may be obtained without any assumptions on 
smoothness of the background fields. The operator pp^ shares some simi- 
larities with the kinetic operator for non-minimally coupled scalars in two 
dimensions (7.12). There are good grounds to believe [305] the functional de- 
terminant (9.11) may describe the spherically symmetric part of the Hawking 
fiux in four dimensions. 

Let us consider the dilatonic operator (7.12). We take ^ = $ for sim- 
plicity^^ . We shall be interested in properties of detD under reflection of the 
dilaton, $ — $. The quantity In det -D($) —In det /}(—$) deflnes the dilaton 
shift under the T-duality transformations in string theory. In the string theory 
context this problem was solved in [378] basing on earlier results of [377] . Here 
we follow presentation of the paper [409]. 



The heat kernel coefficients a2{x,pp^) and a2{x,p^p) are total derivatives. 
Therefore, the second term on the right hand side of (2.33) does not contribute. 

Variation of det D with respect to ^' may be expressed through the scale anomaly. 
Therefore, it is not a problem to find exact dependence of det D on ^' by using the 
methods of sec. 7. 
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Let us introduce two first order differential operators 



and associated second order operators 



Li 



D. 



-e d^e 



LA 



(9.12) 



(9.13) 



We restrict ourselves to fiat two-dimensional manifolds. Therefore, position of 
the indices (up or down) plays no role. 

We calculate again variation of the C, function with respect to variation 

of 



-2sTr 



■s-l 



L,D-- 



-4) 



(9.14) 



The operators in the first term here recombine in D^^. The situation with the 
second term is more subtle. Strictly speaking, to treat this term one has to 
perform the Mellin transform and use analytic continuation in s (see footnote 
26 in sec. 7.1). However, the result of such manipulations is almost obvious: 



L.,.Li 



(9.15) 



where 

H^, = L^Lt (9.16) 

is a matrix operator acting on the space V+ of the vector fields which can 
be represented as v,^ — Lj,/ with some scalar function /. The operator (9.16) 
is not of Laplace type as it contains a complicated differential projector on 

Therefore, even though the variation (9.14) looks very similar to (9.7), for 
example, it cannot be used directly to evaluate ('{0, D^) since we cannot even 
guarantee that (9.15) is regular at s = (in fact, it is not regular, see [409]). 

We also have: 



6C{s, D_) = 2sTr [(5$) (dz' - H'') 



(9.17) 



where the operator 

Hfiu — ^nix'^i^u'L^n'L^' (9.18) 

acts on the space V_ of the vector fields of the form = e^^'Lj^,. All epsilon 
tensors cancel after taking the trace in (9.17) and we arrive at the same formula 
as for but with the replacement $ — $. The spaces K|_ and VI are 
orthogonal, and V-^ + V- = V is the space of all vector fields on M (we neglect 
possible zero modes). Therefore, 



The operator 



2s 



({s, 5$, H + H)- C{s, 5$, D+) - C{s, 5$, D_ 

(9.19) 



(9.20) 
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is of Laplace type. We can now act in the standard manner to obtain 

6 [C'(0, D+) - C{s, D_)] = 2 [a2(<5$, H + H) - a2((5$, D+) - 02(5$, D_) . 

(9.21) 

The right hand side can be evaluated by using (4.27). We leave it is an exercise 
to show that the variation (9.21) is zero (alternatively, one may look up in 
[409]). This leads us to the conclusion that 

C'(0, D+) - C'(0, D_) = - In det D+ + In det D_ = 0. (9.22) 

In the one-dimensional case a similar relation may be obtained by meth- 
ods of supersymmetric quantum mechanics. In this simplest case, and in higher 
dimensions if $ depends on one of the coordinates only, and D_ are 
isospectral up to zero modes and (9.22) follows immediately. In two dimen- 
sions, -D+ and D_ are not isospectral in general. Although and D_ have 
coinciding determinants, other spectral functions can be different. 

Extensions of (9.22) to the case of matrix-valued manifolds with bound- 
aries, and the dilaton-Maxwell theory in four dimensions can be found in [409]. 
A generalisation to arbitrary dimension has been achieved recently [235]. Let 
Ap be a p-form field with the field strength F = dAp. Consider the classical 
action 

jr= [ e-2*FA*F, (9.23) 

where ★ is the Hodge duality operator. Such actions appear, for example, 
in extended supergravities [133] and bosonic M-theory [277]. Instead of the 
dilaton, also a tachyon coupling may appear [300]. We are interested in the 
symmetry properties of the effective action under the transformation p 
n — p — 2, $ — > — $. In higher dimensional supergravity theories this is a part 
of the yS-duality transformation. 

It is convenient to define the twisted exterior derivative 

d^ := e-*de* (9.24) 

and the associated twisted coderivative and twisted Laplacian 

5^ := e*5e-* , A$ := (d$ + . (9.25) 

Since (i| = we have an elliptic complex. The restriction of A$ on the space 
of p forms will be denoted by A^. This twisted de Rham complex was intro- 
duced by Witten in the context of Morse theory and supersymmetric quantum 
mechanics [423]. 

Any p-form can be decomposed as the sum of a twisted exact, twisted 
co-exact, and twisted harmonic form: 

Ap = d^Ap_i + 5^Ap+i + 7p , A$7p = . (9.26) 

The projections on the spaces of twisted exact and twisted co-exact forms will 
be denoted by the subscripts || and ± respectively. 
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We assume that the fields Ap = e~'^Ap have a standard Gaussian measure 
and are to be considered as fundamental fields in the path integral. The action 
given in eq. (9.23) is invariant under the gauge transformation which sends Ap 
to Ap + dfj)Ap_i. This means that the forms which are exact have to be 
excluded from the path integral, but that a Jacobian factor corresponding to 
the ghost fields Ap_i has to be included in the path integral measure. Next we 
note that rf^-exact {p — l)-forms do not generate a non-trivial transformation 
of Ap. Hence, such fields must be excluded from the ghost sector. Then we 
have to include "ghosts for ghosts". This goes on until the zero forms have 
been reached. By giving these arguments an exact meaning, one arrives at 
the Faddeev-Popov approach to quantisation of the p-form actions [175,386, 
349,44]. We note that the procedure of [349,44] is valid also in the presence 
of a dilaton interaction if one simply replaces the ordinary derivatives by the 
twisted ones. As a result, we have the following expression for the effective 
action: 

= ^E(-ir+Mndet(A^U) + W;*°P. (9.27) 

^ fc=0 

depends on certain topological characteristics of the manifold (the Betti 
numbers). We shall neglect W^"^ and some other topological contributions in 
what follows. 

By combinatorial arguments, similar to the presented above, one can show 
that [235]: 

m 

5 {Wp{^) - W^„_p_2(-$)) = ^%) ■ (9-28) 

A;=0 

One should distinguish between 5$ (which is a variation of the dilaton) and 
5$ (which is a twisted co-derivative). Therefore, we have related the variation 
of the effective actions with respect to $ to a combination of the heat kernel 
coefficients which is called the supertrace of the twisted de Rham complex. 
A somewhat surprising feature of the supertrace is that it can be calculated 
for any n, with or without boundaries [235]. For example, the volume term 
in (9.28) does not depend on $ and, therefore, is the standard Euler density 
which is given by 

En := (47r) ~ ^w^^^^ ^m^ji ^'^Ri^j_^j^j^...Ri^_^j_^j^j^_^ (9.29) 

(with n — n/2) for n even. E„ = for n odd. 

It is known [213,44] that for $ = the dual theories are quantum equiv- 
alent: 

l^p(0)-l^„_p_2(0) = 0. (9.30) 
By using this equation as an "initial condition" we can integrate the variation 
(9.28) to obtain: 

Wp{^) - iy„_p_2(-$) = (-1)^ / rf"^v^$E„ . (9.31) 
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Further terms A|,) are of some interest in mathematics. 

They can be calculated [246, 234] for A; < n + 2 and any n. 
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10 Conclusions 



Here we present a short guide in this report. In other words, we are 
going to answer the following question. What should one do if one hkes to 
calculate one-loop counterterms, anomalies, an expansion term in the effective 
action, or something else which is defined by the heat kernel expansion? The 
first step is to find the bulk part of the variation of the classical action (1.2) 
and corresponding operator D. Next one has to bring this operator to the 
canonical form (2.2) with the help of (2.1), (2.3) and (2.4). However, one can 
first consult sec. 3. Probably, relevant expressions can be found there. If the 
problem in question contains neither boundaries nor singularities, one can look 
in sec. 4 for an expression for the heat kernel coefficient, or for a reference, 
or for a method. In the case of the boundaries one has to proceed with sec. 
5, in the case of singularities - with sec. 6. Relations between the heat kernel 
coefficients and quantum anomalies can be found in sec. 7. The results going 
beyond the standard heat kernel expansion arc collected in sec. 8. In particular, 
in this section we explain how one can extract leading non-localities from 
the effective action, and what the heat kernel looks like if the background is 
approximately covariantly constant (in this context we also consider invariant 
operators on homogeneous spaces) . Exact results for the effective action which 
can be obtained with the help of the heat kernel expansion are reviewed in 
sec. 9. 

The Casimir energy is one of the most "classical" applications of the 
heat kernel and zeta function technique (cf. [70]). It follows from the locality 
of the heat kernel expansion that the divergences in the Casimir energy are 
given by volume and surface integrals of some local invariants. Therefore, if 
the boundaries are being moved in such a way that the boundary values of 
the background fields remain unchanged, the boundary contributions to the 
divergences also remain constant. This leads to the well known conclusion [79] 
that there are no boundary divergences in the Casimir force (which is defined, 
roughly speaking as a variation of the Casimir energy under infinitesimal trans- 
lations of the boundary). Consequently, one can assign a unique value to the 
Casimir force (see, e.g., [324]). This observation, however does not mean that 
the quantum field theory on a manifold with boundary is finite. In general, 
some surface counterterms are required (at one-loop order they may be read 
off from sec. 5). Moreover, if the background field are non-trivial, the boundary 
divergences will not be constant. Similar arguments created certain scepticism 
towards reliability of the Casimir energy calculation [252]. This point has not 
been settled so far for a "generic" theory. We may add that in supersymmetric 
theories cancellations between divergences in the bosonic and fermionic sec- 
tors appear if the boundary terms are considered together with the volume 
terms [84], therefore, separation of boundary and volume contributions is not 
always natural for that theories. Renormalization of self-interacting theories 
on manifolds with boundaries was considered in [394,351,317,420,54] where 
one can find further references. Some aspects of the relationship between the 
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Casimir energy calculations and the heat kernel coefficients have been clarified 
recently by Fulling [225]. 

Of course, not everything can be found in this report. Several topics are 
very close to the subject of this review, but are not included. 

(1) The heat kernel expansion can be successfully applied to quantum field 
theory at finite temperature [166] . A new interesting development in this 
field is related to the so called non-linear spectral problem [228] (see [229] 
for an overview). 

(2) The heat kernel expansion has interesting applications to integrable mod- 
els and, in particular, to the KdV hierarchies (see [135] for an elementary 
introduction). 

(3) Recently some attention has been attracted [162,171,384,33] to the so 
called N/D or Zaremba problem which appears when one defines Neu- 
mann and Dirichlet boundary conditions on two (intersecting) compo- 
nents of the boundary. It is unclear whether this problem may have ap- 
plications to quantum theory. 

(4) Instead of considering asymptotics of the heat trace Tr(/e~*^) one can 
also consider an asymptotic expansion for individual matrix elements of 
the heat kernel (/i, e~*^/2) (which are called the heat content asymptotics 
since they remind short time asymptotics of the total heat content in a 
manifold with the specific heat /i and the initial temperature distribution 
/2). Such asymptotics do not contain negative powers of t. More details 
can be found in [297]. 

(5) Many results on the heat kernel asymptotics can be extended to higher- 
order differential operators (see, e.g., [238,268]) and to differential oper- 
ators in superspace [101]. 

Although quite a lot is already known about the heat kernel expansion, 
many interesting problems still remain open. There arc many opportunities to 
extend and generalise the results presented in this report. This refers especially 
to the material of sec. 5-9 where one could add new types of the operators, 
boundary conditions, geometries, and singularities. There is a completely new 
field of research related to the heat kernel expansion where very little has been 
done so far. This is an extension to non-commutative geometry. This problem 
is an especially intriguing one since one can expect very unusual properties 
of the spectral functions because of very unusual properties of corresponding 
field theories in the ultra violet asymptotics. 
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